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Acoustic  Wave  Scattering  from 
a  Random  Ocean  Bottom 

Dajun  Tang 
Abstract 

This  thesis  investigates  low  frequency  acovistic  wave  scattering  from  the  ocean 
bottom.  It  is  divided  into  four  parts.  The  first  part  models  the  ocean  bottom  as  a 
fluid  medium  where  sound  speed  and  density  are  constants  except  a  layer  in  which 
the  density  is  still  a  constant,  but  the  sound  speed  is  composed  of  a  large  constant 
and  superimposed  with  a  small  random  component.  It  is  assumed  that  the  random 
soimd  speed  is  horizontally  well  correlated  and  vertically  poorly  correlated. 

In  the  second  part,  an  integro-differential  equation  method  is  used  to  calculate 
the  scattering  from  that  random  layer.  Emphasis  is  put  on  the  study  of  the  spatial 
correlation  of  the  scattered  field.  It  is  found  that  the  spatial  correlation  length  of 
the  scattered  field  is  related  to  the  correlation  length  of  the  scatterer,  therefore 
it  is  possible  to  invert  the  bottom  correlation  length  by  measuring  the  spatial 
correlation  of  the  scattered  field  using  multiple  receivers.  Also  included  in  this 
part  are  an  estimation  of  energy  loss  in  the  coherent  field,  a  discussion  on  the 
influence  of  bottom  anisotropy,  and  a  comparison  between  the  integro-dilferential 
equation  method  and  the  Bom  approximation. 

The  third  part  concerns  the  influence  of  the  bottom  roughness.  A  small  rough¬ 
ness  is  added  to  the  water/bottom  interface  and  perturbation  method  is  used  to 
calculate  the  scattering  from  the  roughness.  Under  the  assumption  that  the  rough¬ 
ness  and  the  volume  inhomogeneity  are  uncorrelated,  compau-isons  are  made  on 
the  scattering  strengths  between  roughness  scattering  and  volume  scattering,  and 
the  spatial  correlation  of  the  total  scattered  field  is  evaluated.  It  is  found  that  at 
low  frequencies,  volume  scattering  cannot  be  ignored  when  the  seafloor  is  not  very 
rough,  and  it  is  possible  to  distinguish  the  two  scattering  mechanisms  by  measuring 
spatial  correlations  of  the  scattered  field. 

In  the  fourth  part,  the  classical  definition  of  the  bottom  scattering  cross  section 
is  questioned,  and  a  new  set  of  parameters,  the  Scattering  Correlation  Coefficient, 
is  introduced,  which  does  not  have  the  ambiguities  associated  with  the  conventional 
notion  of  the  scattering  cross  section. 
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Chapter  1 


INTRODUCTION 


Sound  scattering  Is  a  very  important  part  of  tinderwater  acoustics  research. 
Since  World  War  11,  considerable  effort  has  been  put  into  understanding  the  mech¬ 
anisms  that  cause  scattering  of  sound  waves  which  propagate  in  the  ocean[l].  The 
ocean,  as  an  acoustical  system,  is  extremely  complicated.  Scattering  can  happen 
almost  everywhere  in  the  ocean  and  it  covers  a  very  wide  bzmd  of  frequencies  (from 
seismic  waves  of  below  one  Hz  up  to  MHz). 

Conventionally,  acoustic  scattering  can  be  categorized  into: 

a)  Ocean  surface  scattering:  here,  scattering  is  caused  by  the  rough  ocean 
surface  and  a  thin  bubble  layer  immediately  below  the  surface. 

b)  Volume  scattering  within  the  water:  here,  scattering  is  due  to  organic  species 
or  temperature  and  density  fluctuations. 

c)  Bottom  scattering:  this  part  of  scattering  is  from  the  water /bottom  rough 
interface  and  the  inhomogeneities  (fluctuations  of  sound  speed  and  density,  isolated 
rocks,  etc.)  within  the  bottom. 

Each  of  these  categories  is  further  divided  into  sub>categories  according  to  the 
size,  nature  and/or  the  petition  of  the  scatterers,  or  the  sound  frequency  being 
used.  In  this  research,  all  the  efforts  will  be  concentrated  on  the  bottom  scattering 
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part,  although  some  of  the  methods  and  results  are  applicable  to  volume  scattering 
in  the  water  column  as  well. 


1.1  Motivation 

Ocean  bottom  scattering  is  a  process  of  considerable  importance  in  the  oper¬ 
ation  of  underwater  acoxistic  systems.  In  some  cases,  bottom  scattering  is  a  major 
source  of  interfering  reverberation  that  restricts  sonar  performance;  in  other  cases, 
it  provides  an  effective  means  of  remotely  measuring  the  physical  properties  of  the 
ocean  bottom.  For  basic  understanding  of  the  creation  and  distribution  of  the 
bottom  scatterers,  please  see  Kennett[2]. 

A  great  deal  of  research  has  been  concentrated  on  water /bottom  rough  inter¬ 
face  scattering.  A  number  of  bottom  roughness  models  have  been  developed  over 
the  years  to  estimate  the  magnitude  of  either  backscatter  or  forwardscatter  from 
the  seabed  boundary,  as  a  ftinction  of  incident/scattering  angles,  frequency,  bottom 
type  and  roughness[3,4,5,6,7,8].  All  these  models  are  developed  by  modifying  rough 
ocean  surface  scattering  theories  which  impose  a  pressure-release  boundary  con¬ 
dition.  Their  success  is  limited  because  the  scatterers  beneath  the  water /bottom 
interface  are  ignored.  Especially  when  the  bottom  is  smooth,  the  roughness  models 
fail  to  predict  the  correct  scattering  strength  and  frequency  dependence[9, 10,13]. 

It  is  generally  accepted  that  volume  scattering  caused  by  the  inhomogeneities 
beneath  the  water/bottom  interface  cannot  be  heglected[ll,12,9,13].  In  some  cases, 
when  the  interface  is  smooth,  volume  inhomogeneities  are  the  dominant  contributor 
to  scattering  [8,9,13].  A  number  of  researchers  have  developed  voltune  scattering 
theories  for  the  bottom.  Stockhausen|14]  developed  a  scattering  model  for  a  bot¬ 
tom  consisting  of  a  uniform  set  of  spherical  particles,  in  which  the  scatterers  are 
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not  correlated.  Ivakin  and  Lysanov  [10]  and  Ivakin[15]  employed  Chemov’s[16] 
work  on  wave  propagation  in  inhomogeneous  media  to  develop  bottom  scatter¬ 
ing  models.  They  do  not  include  the  lateral  wave  contributions.  Hines[9,13]  fol¬ 
lowed  a  similar  approach  but  included  lateral  wave  contributions.  He  compared  his 
model  with  his  own  laboratory  measurements  and  with  other  available  data.  Both 
Ivakin/Lysanov’s  model  2md  Hines*  model  considered  the  scatterers  in  the  bottom 
as  correlated  random  quantities,  and  this  approach  appears  to  be  appropriate  after 
comparing  it  with  high-frequency  experimental  data. 

Efforts  have  also  been  made  to  incorporate  roughness  and  volume  inhomo¬ 
geneities  into  a  single  model.  Ivakin  and  Lysanov[17]  include  refraction  at  a  ran¬ 
domly  rough  interface  in  their  sediment  scattering  model.  Refraction  by  large-scale 
roughness  is  considered,  but  diffraction  by  the  smaller  scales  is  neglected.  Jackson 
et  a/. [12]  developed  a  composite  roughness  model  for  high-frequency  backscatter- 
ing,  where  the  volume  scattering  is  converted  into  equivalent  surface  scattering 
under  the  assumption  that  high  absorption  rates  prevent  acoustic  energy  from 
penetrating  deep  into  the  bottom. 

However,  problems  in  bottom  scattering  still  exist.  First,  none  of  the  existing 
volume  scattering  theories  considers  the  coherent  sound  energy  loss  due  to  scatter¬ 
ing  while  propagating  in  the  random  medium,  and  only  absorption  of  sound  due 
to  heat  conversion  is  included.  When  the  acoustic  frequency  is  so  high  that  the 
heat  absorption  is  much  greater  than  the  scattering  loss,  these  theories  work  well. 
For  low  frequencies,  however,  ignoring  volume  scattering  loss  is  not  justified.  This 
effect  is  yet  to  be  addressed  in  the  literature. 

Secondly,  when  taking  bottom  inhomogeneities  as  a  random  process,  the  under¬ 
standing  of  the  auto-correlation  function  of  this  process  is  essential  for  modeling 
this  process.  The  existing  theories  concentrate  only  on  finding  the  acoustic  scat¬ 
tering  strength,  without  going  further  to  examine  the  spatial  correlation  of  the 
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scattered  acoustic  field.  The  acoxistic  scattering  strength  (standard  deviation  of 
the  randomly  scattered  field)  does  not  give  a  complete  characterization  of  the  sta¬ 
tistical  properties  of  the  wave  field.  The  statistical  properties  of  the  fiuctuations 
of  the  wave  field  can  be  characterized  more  completely  by  using  a  correlation  func¬ 
tion  of  the  scattered  field.  However^  the  question  of  how  to  find  and  utilize  the 
correlation  of  the  scattered  field  in  the  bottom  scattering  process,  which  is  of  both 
theoretical  and  practical  interest,  has  not  yet  received  much  attention. 

Thirdly,  it  has  been  argued  that  when  the  water/bottom  interf3u:e  is  not  very 
rough,  the  voliime  scattering  contribution  to  the  total  bottom  scattering  process 
will  become  important.  However,  there  is  no  quantitative  description  of  this  effect. 
More  precisely,  given  the  ramdom  parameters  of  the  volume  inhomogeneity,  what 
degree  of  interface  roughness  will  result  in  an  equivalent  aimount  of  scattering  for 
a  given  frequency  range? 

Fourthly,  bottom  anisotropy  (i.e.,  the  bottom  random  parameters  have  different 
values  if  measured  along  different  directions)  is  yet  to  be  investigated.  When 
the  bottom  random  processes  is  not  isotropic,  what  impact  it  will  have  upon  the 
scattered  field? 

Finally,  all  existing  volume  scattering  models  use  a  far-field  assumption,  and 
the  scattering  cross  section  or  the  scattering  coefficient  is  used  to  characterize  the 
scatterers.  However,  for  low-frequency  scattering  in  a  layered  random  medi\im, 
both  near-  and  far-field  solutions  are  necessary,  and  the  classical  definition  of  scat¬ 
tering  cross  section  is  no  longer  appropriate. 

So  far,  almost  all  bottom  scattering  experiments  are  conducted  by  using  a  sin¬ 
gle  soTirce  and  a  single  receiver.  A  typical  example  can  be  found  in  Merklinger’s 
paper[ll].  In  his  experiment,  the  source  and  the  receiver  are  placed  along  a  vertical 
line  and  are  arranged  as  close  together  as  possible  to  ensure  monostatic  backscat- 
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tering  geometry.  By  assuming  the  only  scatterers  axe  along  the  water/bottom 
interface,  and  that  these  scatterers  are  homogeneously  distributed,  he  finds  the  in> 
sonified  area,  an  annulus  controlled  by  the  source  duration  and  the  sound  speed  in 
the  water,  and  calculates  the  backscattering  strength  as  a  function  of  incident  an¬ 
gle.  He  noticed  that  there  are  some  features  in  the  data  which  cannot  be  explained 
by  interface  roughness,  and  he  speculates  that  these  featiires  are  caused  by  sub¬ 
bottom  inhomogeneities.  However,  single-source-single-receiver  experiments  of  this 
type  cannot  distinguish  interface  scattering  from  sub-bottom  volume  scattering. 

Because  of  the  fast  development  of  technology,  it  is  now  feasible  to  use  acoustic 
arrays  to  measure  bottom  scattering.  Frisk[l8]  proposes  to  use  vertical  arrays  to 
measure  low-frequency  bottom  scattering.  Motivated  by  Frisk’s  proposail,  this  the¬ 
sis  intends  to  address  the  aforementioned  problems.  An  integro-differential  equa¬ 
tion  method  will  be  used  to  formulate  the  problem.  The  major  contributions  of 
this  thesis  are: 

1  The  coherent  energy  loss  due  to  scattering  is  estimated.  As  a  result,  the 
scattering  theory  developed  in  this  thesis  improves  upon  the  Bom  approximation. 

2  For  the  first  time  in  bottom  scattering  problems,  spatial  correlation  of  the 
scattered  field  is  studied  in  detail  for  two  fluid  half-spaces  with  the  lower  half-space 
being  random.  This  extends  our  imderstanding  of  the  bottom  scattering  process 
and  a  technique  to  invert  random  parameters  in  the  bottom  is  proposed. 

3  Scattering  from  a  rough  interface  on  top  of  a  layer  of  volume  inhomogeneity 
is  studied,  amd  a  compairison  between  rough  interface  scattering  and  scattering  due 
to  volume  inhomogeneities  is  made.  We  handled  the  interface  scattering  pau't  by 
using  the  boundary  perturbation  theory[19]  and  assumed  the  roughness  and  the 
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volume  inhomogeneity  axe  uncorrelated.  The  result  of  this  analysis  suggests  the 
possibility  of  distinguishing  the  two  scattering  mechanisms  experimentally. 

4  The  methods  used  in  this  thesis  allow  the  inclusion  of  anisotropy  of  the  scat- 
terers,  and  so  the  anisotropic  influences  on  scattering  are  also  discussed  at  length. 

5  The  Scattering  Correlation  Coefficient  is  introduced,  which  clarifies  the  am> 
biguities  of  using  the  concept  of  scattering  cross  section  in  bottom  scattering  prob¬ 
lems.  A  formulation  is  found  based  on  the  Scattering  Correlation  Coefficient  to 
calculate  the  scattering  field  due  to  a  point  source.  It  will  help  tis  understand 
conceptually  the  physical  process  of  volume  scattering  in  complicated  media,  and 
it  can  be  used  as  the  starting  point  to  solve  the  near-field  problem. 


1.2  Review  of  Relevant  Literature 

The  purpose  of  this  section  is  not  to  provide  an  exhaustive  review  of  all  the 
existing  literature  on  bottom  scattering.  Since  this  research  involves  several  differ¬ 
ent  areas,  it  is  intended  only  to  discxiss  the  references  which  are  directly  relevant 
to  each  area. 

1.2.1  On  wave  propagation  in  a  random  medium 

There  are  three  excellent  books,  one  by  Chemov[16],  one  by  Tatarskii[20],  and 
the  third  by  Ishimaru[22],  which  contain  systematic  treatments  of  the  theory  of 
wave  propagation  in  a  medium  with  random  inhomogeneities.  They  include  the 
study  of  wave  propagation  xising  the  ray  approximation;  the  diffraction  theory  of 
wave  propagation  and  related  approximation  methods;  the  application  of  quantum 


15 


field  theory;  the  infiuences  of  fiuctuations  in  the  incident  wave  on  the  diffraction 
image  formed  by  a  focusing  system;  and  some  experimental  data  analyses.  Many 
basic  concepts  in  propagation  and  scattering  in  a  random  medium  are  clearly  de¬ 
scribed.  Since  the  authors  treat  the  problem  of  wave  propagation  as  a  general 
physical  phenomenon,  the  books  cover  almost  all  branches  of  wave  motion.  How¬ 
ever,  the  problem  of  wave  propagation  in  a  random  medium  with  discontinuous 
interfaces  is  not  included  since  this  is  still  an  on  going  research  subject. 

Also,  there  is  a  review  paper  by  Barabanenkov  tt  al.  [40]  on  the  ‘Status  of  the 
theory  of  propagation  of  waves  in  a  randomly  inhomogeneous  medium” ,  where  the 
authors  restrict  their  review  to  volume  scattering  in  media  without  discontinuous 
interfaces.  The  method  of  small  perturbations,  and  scattering  by  large-scale  inho¬ 
mogeneities  are  discussed  in  detail.  There  is  also  a  very  insightful  introduction  to 
the  theory  of  multiple  scattering.  Literature  cited  ranges  from  scattering  in  radio 
physics,  optics,  acoustics  to  scattering  in  plasma  physics  and  other  branches  of 
physics. 

Flatt4[23,24]  and  his  colleagues  have  performed  extensive  studies  on  acoiistic 
waves  propagating  in  a  temporally  and  spatially  varying  oceein,  with  emphasis  on 
the  effect  of  internal  waves. 

1.2.2  On  volume  scattering 

The  Born  approximation  is  the  most  widely  used  method  in  underwater  acous¬ 
tics  in  solving  weak  volume  scattering  problems.  The  original  form  of  the  Born 
approximation  [64]  is  for  those  scattering  problems  that  the  medium  is  homoge¬ 
neous  except  for  the  scatterers.  Later  the  Bom  approximation  is  extended  to  cover 
the  case  that  the  background  medium  is  stratified  [  10).  We  will  call  it  the  Extended 
Born  approximation  in  this  thesis.  When  the  background  medium  is  stratified,  both 
the  background  wave  and  the  scattered  wave  become  more  complicated  because  of 
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multi-path  effect. 

Ivakin  and  Lysanov[lO,25,15]  followed  Chernov’s  theory[l6],  assuming  that 
the  random  scattering  is  caused  by  ’‘sharply  anisotropic  random  inhomogeneities  ( 
fluctuations  of  the  refractive  index):  large-scale  in  the  horizontal  plane  and  small- 
scale  in  depth”.  They  used  the  Extended  Bom  approximation  and  the  far-field 
assumption  (steepest-descent  method)  to  obtain  a  backscattering  formula.  The 
authors  explained  the  distinctive  featxires  observed  in  the  frequency  and  angular 
dependence  of  the  backscattering  coefficient  at  small  grazing  angles.  Although 
the  water/bottom  interface  discontinuity  is  accounted  for  by  using  the  Extended 
Born  approximation,  the  authors  ignored  the  interface  wave  contribution  in  their 
calculations. 

Hines’  Ph.D.  thesis[9]  in  1988  is  a  systematic  study  of  high-frequency  acoustic 
wave  backscattering  from  volume  inhomogeneities.  It  contains  an  extensive  lit¬ 
erature  survey  on  wave  scattering  both  from  rough  interfaces  and  from  volume 
inhomogeneities.  Hines’  theory  also  followed  Chernov’s  approach  [16],  and  the  Ex¬ 
tended  Born  approximation  and  the  far-field  approximation  is  used  to  calculate 
volume  backscattering.  His  assumption  of  bottom  randomness  is  that  both  the 
soimd  speed  and  the  density  fluctuate,  and  the  two  are  perfectly  correlated  through 
medium  porosity.  The  correlation  function  he  \ised  is  exponential.  One  step  fur¬ 
ther  than  Ivakin  and  Lysanov,  Hines  also  included  the  interface  wave  contribution. 
The  results  are  compared  to  his  own  high  frequency  laboratory  measurements  as 
well  as  published  experimental  data  conducted  by  other  researchers. 

1.2.3  On  the  integro-differential  equation  method 

This  method  is  the  so-called  “bi-local”  approximation  to  the  Dyson  equation 
for  calculating  the  coherent  part  of  the  acoustic  field.  (For  reference  see  Tatarskii, 
Bourret  and  Keller  [20,26,27]).  This  approach  broadens  the  domain  of  validity  of 
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the  single-scattering  (Born)  approximation  to  a  considerable  extent.  Most  of  the 
research  is  for  the  case  of  an  unbounded  random  medium.  When  a  random  medium 
is  confined  to  a  finite  region,  it  is  a  much  more  difficult  mathematical  problem. 
Rosenbaum[28]  studied  the  interface  effects  and  concluded  that  the  introduction 
of  interfaces  will  change  the  solution  of  the  integro-differential  equation.  For  a 
one-dimensional  random  medium,  the  interface  perturbs  the  coherent  wave  motion 
to  a  significant  extent  over  a  large  distance  from  it;  for  the  three  dimensional 
model  which  Rosenbaum  used,  the  interface  effects  are  confined  to  a  well-defined 
'transition’'  layer  of  the  order  of  a  wavelength.  However  he  only  considered  the  case 
of  normal  incidence.  Kupiec  et  a/. [29]  studied  the  reflection  for  a  special  medium 
where  the  average  refractive  index  in  the  random  portion  is  equal  to  the  refractive 
index  in  the  homogeneous  region. 

Kurtepov[30]  demonstrated  one  of  the  very  few  applications  of  the  integro- 
differential  equation  method  in  underwater  acoustics.  He  studied  the  influence  of 
internal  waves  on  shallow  water  propagation.  Assuming  that  both  bovmdaries  of  the 
shallow  water  waveguide  are  impenetrable,  he  obtained  the  attenuation  coefficients 
for  each  propagating  normal  mode. 

1.2.4  On  the  rough  interface  model 

Due  to  the  volume  of  work  done  on  this  subject,  it  b  virtually  impossible  to 
outline  all  related  papers.  A  good  review  paper  by  Fortuin[31]  summarizes  pressure- 
release  rough  surface  scattering  theories  up  to  1969.  When  the  roughness  height 
is  relatively  small,  compared  with  the  acoustic  wavelength,  perturbation  methods 
are  widely  used.  Kuperman[6]  employed  the  perturbation  method  developed  by 
Bass[32]  for  the  rough  sea  surface,  and  obtained  a  set  of  boundary  conditions 
to  calculate  the  reflection  and  transmission  at  a  two-fluid  rough  interface.  Later 
Kuperman  and  Schmidt[7]  extended  the  model  to  handle  elastic  wave  scattering 
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in  a  horizontally  stratified  ocean.  Theories  on  large  amplitude  roughness  will  not 
be  cited  here  since  this  thesis  studies  only  the  case  where  the  interface  is  relatively 
smooth.  Recently,  there  are  researchers  who  model  the  rough  ocean  bottom  as 
fractals  [33,34,35,36,37].  One  major  advantage  of  this  approach  is  that  if  the  rough 
boimdary  is  indeed  fractal,  the  self-similar  property  of  the  fractal  will  enable  us  to 
just  measure  the  statistics  of  the  roughness  at  one  scale  and  predict  the  statistics  of 
the  roughness  on  all  other  scales.  In  this  thesis,  Bass  and  Fuks’[l9]  single  scattering 
model  will  be  used  as  the  means  to  compare  the  relative  scattering  strengths  of 
volume  and  rough  interface  scatterings. 

1.2.5  On  spatial  coherence  of  scattered  field 

In  underwater  acoustics,  there  has  been  of  little  study  on  the  spatial  coherence 
of  the  bottom  scattered  field.  However,  research  on  spatial  coherence  of  rough 
surface  scattered  fields  exists.  Notably,  Clay  and  Medwin  [42,43]  in  1970  studied 
the  dependence  of  the  spatial  and  temporal  correlation  of  forward-scattered  sound 
on  the  rough  surface  statistics.  They  used  the  Kirchhoff  approximation  to  formulate 
the  problem,  and  related  the  spatial  and  temporal  correlations  of  the  acoustical 
field  to  the  statical  description  of  the  rough  surface.  In  a  laboratory  experiment, 
a  wind-driven  sea  surface  was  simulated,  and  the  spatial  and  temporal  correlation 
of  the  acoustic  field  was  measured  and  compared  with  their  theory. 

Kinney  and  Clay[44]  studied  spatial  coherence  by  using  the  facet-ensemble 
method  in  1984.  They  compared  the  experimental  result  by  Medwin  and  Clay[43] 
with  Eckart  theory  and  the  facet-ensemble  method,  and  found  good  agreement  be¬ 
tween  the  experiment  data  and  the  theory.  In  another  paper[45],  the  same  authors 
studied  if  the  spatial  wavenumber  power  spectrum  of  a  rough  surface  is  sufficient 
to  predict  the  amplitude  and  spatial  coherence  of  acoustic  energy  scattered  from 
that  surface,  and  the  conclusion  was  negative. 
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1.3  Outline  of  the  Content  of  Each  Chapter 

Chapter  2  defines  the  volume  scattering  scenario.  Assuming  the  scattering  is 
ca\ised  by  the  random  variation  of  sound  speed  in  a  layer  beneath  the  water/bottom 
interface,  the  wave  equation  is  split  into  two  coupled  equations,  governing  the 
coherent  component  and  the  random  component  of  the  total  field.  The  relationship 
between  these  two  components  is  discxissed. 

Chapter  3  is  devoted  to  solving  for  the  coherent  wave  field.  An  incident  plane 
wave  is  assumed,  and  the  coherent  refiection  coefficient  between  two  fiuid  half¬ 
spaces  is  obtained  for  all  incident  angles  and  different  frequencies.  Different  back¬ 
ground  parameters  will  be  xised  in  the  calculation,  and  the  results  are  compau’ed 
to  the  reflection  coefficients  when  no  volume  inhomogeneities  are  present,  but  heat 
absorption  is  included. 

Chapter  4  calculates  the  randomly  scattered  field,  and  studies  the  spatial  cor¬ 
relations  (vertical  and  horizontal)  of  this  field.  It  is  found  that  the  spatial  correla¬ 
tion  of  the  scattered  field  depends  strongly  on  the  correlation  length  of  the  bottom 
random  sound  speed  variation;  as  a  result,  by  sweeping  frequencies,  the  bottom 
correlation  length  can  be  determined  experimentally.  The  results  are  compared 
with  that  of  the  Extended  Born  approximation,  and  some  conclusions  are  made 
as  to  the  applicability  of  the  Extended  Bom  approximation.  Finally,  the  influence 
of  bottom  anisotropy  is  examined.  If  the  bottom  correlation  function  has  a  longer 
correlation  length  in  one  direction  and  a  shorter  one  in  the  orthogonal  direction,  it 
is  found  that  the  longer  and  the  shorter  correlation  length  of  the  bottom  correla¬ 
tion  function  are  both  sensitive  to  the  acoustic  frequency  and  the  orientation  the 
the  receiving  arrays. 

Chapter  5  compares  the  differences  between  scattering  from  volume  inhomo¬ 
geneities  emd  scattering  from  small  interface  roughness.  We  consider  for  various 
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parameters,  which  of  the  two  mechanisms  will  be  dominant.  When  either  is  dom¬ 
inant,  the  spatial  correlation  of  the  combined  scattered  field  is  investigated,  and 
the  possibility  of  using  it  to  distinguish  the  two  mechanisms  is  discussed. 

In  Chapter  6,  the  classical  definition  of  scattering  cross  section  of  ocean  bottom 
scattering  is  questioned.  A  new  set  of  parameters,  the  "Scattering  Correlation 
Coefficients”,  is  defined  which  does  not  have  the  ambiguities  associated  with  the 
conventional  notion  of  scattering  cross  section.  The  mathematical  formulation  for 
point  source  scattering  is  given  in  detiul.  The  significance  of  the  formulation  and 
the  difiiculties  in  solving  this  problem  are  examined. 

Chapter  7  concludes  the  thesis  and  addresses  future  research  works  related  to 
this  thesis.  What  are  the  necessary  geological  parameters  to  be  used  as  inputs  in 
the  theory  developed  in  the  thesis  are  discussed  in  detail. 
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Chapter  2 


FORMULATION 


2.1  The  Integro-differential  Equation 

In  the  thesis,  we  are  interested  in  acoustical  wave  scattering  by  an  inhomoge¬ 
neous  ocean  bottom.  The  medium  under  study  is  schematically  shown  in  Fig.2.1. 

The  upper  half-space,  z  >  0,  is  composed  of  homogeneous  water  having  con¬ 
stant  density  p\  and  constant  sound  speed  Cx.  The  lower  half-space,  z  <  0,  is  a 
homogeneous  fluid  ocean  bottom  having  constant  density  ps  and  constant  sound 
speed  cj,  except  for  the  layer  between  z  =  0  and  z  =  — h,  where  the  medium  has 
the  same  constant  density  pj  as  the  rest  of  the  bottom,  but  the  sound  speed  is 
composed  of  a  constant  C]  and  superimposed  with  a  random  component  Sc,  where 
^c  is  a  fimction  of  position  (x,y,z): 

Cuytr  =  Cj  +  ^c(l,y,  z).  (2.l) 

We  also  assume  {Sc[x,y,x))  =  0,  and  so  (e<av*r)  =  cs,  where  the  brackets  (•  •  •) 
indicate  the  ensemble  average.  The  thickness  of  the  random  layer  is  a  free  param¬ 
eter,  which  can  be  chosen  as  large  as  infinity.  The  coherent  reflection  and  random 
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scattering  from  this  random  layer  is  the  primary  concern  in  this  thesis. 

The  validity  of  the  theory  that  will  be  developed  is  not  restricted  to  the  layered 
fluid  media  as  shown  in  Fig.2.1.  The  bottom  can  be  any  fluid,  fluid-elastic,  or 
elastic  layered  media,  and  the  random  layer  is  not  necessarily  immediately  beneath 
the  water;  there  can  be  other  homogeneous  layers  beneath  the  water  and  above 
the  random  layer.  Further,  more  than  one  random  layer  b  allowed.  For  a  more 
complete  model,  attenuation  of  both  compressional  waves  and  shear  waves  should 
be  included.  The  attenuation  can  be  incorporated  into  the  present  theory  without 
difficulty.  Another  important  issue  is  to  model  the  background  sound  speed  proflle 
in  the  random  layer  with  a  gradient,  which  is  closer  to  a  real  sediment  sound  speed 
proflle[63].  This  profile  will  allow  the  incident  wave  to  refract  back  into  the  water 
column,  hence  causing  additional  scattering. 

However,  in  this  thesis  we  will  restrict  the  study  to  the  case  shown  in  Fig.2.1, 
emphasizing  the  basic  physics  of  scattering  from  a  simple  random  layer.  This  will 
allow  us  to  closely  observe  the  consequences  of  the  random  layer  on  the  behavior 
of  the  acoustic  wave  scattering. 

Generally,  the  density  of  the  random  layer  should  be  treated  as  a  random  pa¬ 
rameter  as  well.  We  must  therefore  know  the  correlation  between  the  random 
sound  speed  and  the  random  density.  However,  there  is  little  experimental  data 
available  to  extract  such  a  relationship.  We  will  only  study  the  scattering  by  sound 
speed  fluctuations  in  this  thesis.  The  effects  of  density  fluctuations  can  be  incor¬ 
porated  into  the  present  theory  when  the  sound  speed-density  correlation  is  known. 

This  theory  will  be  suitable  for  studying  scattering  by  a  sedimented  bottom, 
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where  the  interface  roughness  does  not  necessarily  dominate  the  scattering  process, 
so  volume  scattering  can  be  important;  also  the  inhomogeneities  can  be  treated  as 
small  perturbations  since  the  scatterers  are  due  to  inhomogeneous  sedimentation, 
where  the  sound  speed  variation  is  not  very  large.  Scattering  by  rocky  bottoms 
and  bottoms  with  ocean  ridges  should  be  treated  differently. 

For  the  remainder  of  this  thesis,  we  shall  use  the  following  subscript  notations. 
We  use  1,  2,  and  3  to  respectively  designate:  the  water  (z  >  0),  the  random 
layer  (0  >  z  >  — h),  and  the  homogeneous  half-space  beneath  the  random  layer 
(z  <  —h).  With  suppressed  time  dependence  e“*"*,  the  pressure  wave  equation  for 
the  homogeneous  water  medium  is: 


V’Pi  +  klPi  =  0.  (2  >  0) 


(2.2) 


Likewise,  the  pressure  wave  equation  for  the  homogeneous  bottom  medium  {z  < 
—h)  is: 


V’P,  -h  kiPj  =  0,  {z<  -h) 


(2.3) 


where  fci  =  w/ci  is  the  wavenumber  in  the  water,  and  fcj  =  u/c2  is  the  wavenumber 
in  the  homogeneous  bottom.  We  note  that  the  average  wavenumber  in  the  random 
layer  equals  the  wavenumber  of  the  homogeneous  bottom.  The  pressure  wave 
equation  in  the  random  layer  is: 


V’p,  +  (- 


yp,  =  0.  (0  >  z  >  -h) 


(2.4) 


Cj  +  Sc^ 

We  assume  the  random  part  of  the  sound  speed  is  small,  6c/ C2  <.  1,  the  wavenum¬ 
ber  in  the  random  layer  will  take  on  the  form: 

UJ 


^  —d  _ 

Cj  -h  Sc  Cj  Cj 


Hence 


,  w  w*.  26c.  ,  ,, 


(2.5) 


(2.6) 
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where  €{x,y^z)  =  — 25c/ca. 

Therefore,  Eq.(  2.4)  becomes: 

V’Pj  +  fc|[l  +  €(z,y,z)]Pj  =  0.  (0  >  z  >  -h)  (2.7) 

Since  there  is  no  general  solution  to  the  above  Helmholtz  equation  with  statistically 
random  coefficient  e,  we  decompose  Pi  into  a  coherent  pairt  {Pi)  and  a  scattered 
part  Pi  throughout  the  media: 

Pi^{P^-\-pl  (i  =  1,2,3)  (2.8) 

The  wave  equations  governing  (Pi),  pj,  (Pj),  and  pj  are  the  same  homogeneous 
Helmholtz  equations  as  Eqs.(  2.2)  and(  2.3),  whereas  the  wave  equation  for  the 
random  layer,  Eq.(  2.7),  becomes  a  pair  of  coupled  differential  equations  governing 
the  coherent  part  (Pj)  and  the  scattered  part  pj  respectively: 


V'(P,)  +  k\{P,)  =  -*|{ep5), 

(2.9) 

VV5  +  *jPi  =  -^k\{Pt)  +  tlKvS)  -  epj). 

(2.10) 

These  are  inhomogeneous  Helmholtz  equations  with  unknown  source  terms  to  the 
right  of  the  equalities.  Notice  that  the  term  in  square  brackets  of  Eq.(  2.10)  is  of 
second-order  smallness.  By  eliminating  that  term  (called  bi-local  approximation[20]), 
we  see  that  the  coherent  field  determines  the  source  strength  of  the  scattered  field 
in  Eq.(  2.10);  in  turn,  the  scattered  field  determines  the  source  strength  of  the 
coherent  field  in  Eq.(  2.9).  Applying  Green’s  Theorem  to  Eq.(  2.10),  we  have[20] 

p;(fi)  =  ^//.,/ (211) 

where  v*  is  the  random  region,  and  R  and  R!  indicate  source  position  (z',y',  s') 
and  receiver  position  (z,y,z),  respectively.  G  is  the  Green’s  function  which  will 
be  given  in  detail  later.  Substituting  Eq.(  2.11)  into  Eq.(  2.9),  we  find  that  the 
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coherent  part  of  the  pressure  field  in  the  random  layer  fulfills  the  following  integro- 
difierential  equation: 

-  (*}/4t)  f  \^)dv'.  (2.12) 

The  above  equation,  along  with  the  coherent  fields  in  media  1  and  3  and  appropriate 
boimdaiy  conditions,  determines  the  coherent  field  in  the  random  layer.  After 
finding  (Pi),  we  will  be  able  to  use  £q.(  2.11)  to  calculate  the  scattered  field. 
Notice  that  in  £q.(  2.11),  if  we  substitute  Pf,  the  pressure  field  in  the  absence  of 
6,  instead  of  using  (Pi),  then  it  becomes  the  Extended  Bom  approximation[9]. 

£q.(  2.12)  has  been  studied  for  cases  without  velocity  discontinuities  by  a  num¬ 
ber  of  authors  for  various  kinds  of  isotropic  correlation  functions[20,28,29].  In  the 
present  problem,  there  is  a  sound  speed  discontinuity  at  the  water /bottom  inter¬ 
face.  The  presence  of  this  discontinuity  will  result  in  a  more  complicated  Green’s 
Function  compared  with  that  when  the  discontinuity  is  absent.  The  way  the  inter¬ 
face  problem  is  handled  will  be  given  in  Chapter  3.  In  the  rest  of  this  chapter,  we 
will  introduce  the  Green’s  Function  G  and  the  correlation  function  (£(M)e(B')}. 


2.2  The  Green’s  Function 

Since  we  have  chosen  that  {c^ytr)  —  ci,  the  backgroimd  sound  speed  is  just  two 
homogeneous  half-spaces.  The  Green’s  Function  for  two  half-spaces  with  source 
located  in  the  lower  half-space  satisfies: 


(V»  +  k\)Gi(R  1  5^)  =  0, 

2  >  0 

(2.13) 

(V’  +  *|)G,(5  I  R')  =  -4)ri(S  -  R% 

2  <  0 

(2.14) 
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and  the  boundary  conditions  at  z  0  are: 


I  I'.y',*')  =  Gj(i,»,0  I 


(*•15) 


(2.16) 


The  Green’s  function  for  the  lower  half-space  can  be  expressed  in  the  following 
integral  form[54]: 

G2[R  1  B')  =  y*  y  'I  -I-  (2.17) 

where  vsi  is  the  plane  wave  reflection  coefficient  with  the  plane  wave  incident  from 
below  the  water/bottom  interface: 


Vjl  —  i  i«.a.oi 

and  and  ^2  are  the  vertical  wave  numbers  in  the  upper  and  lower  half-spares 
respectively: 


(2.18) 


6  = 
6  = 


(2.19) 

(2.20) 


and  /  is  the  vector  connecting  the  horizontal  coordinates  of  the  source  and  the 
observation  points: 


r=  r  -  r"  =  (x  -  x')r-|-  (y  -  y')j. 


(2.21) 


Notice  that  the  Green’s  function  for  a  medium  without  any  discontinuity  would  be 
the  same  as  Eq.(  2.17)  except  that  U2x  =  0.  The  inclusion  of  t>2i  in  the  Green’s 
function  makes  the  solution  of  the  coherent  fleld  much  more  difficult  to  find  than 
in  the  case  when  V21  =  0. 
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2.3  The  Correlation  Function 


We  assume  that  the  correlation  ftmction  of  the  quantity  e(J2)  has  the  following 
stationary  and  anisotropic  form,  wMch  is  horizontally  and  vertically  decoupled[10,50]: 

<c(B)e(B'))  =  =  4<T>V>r(i)Af  (|r  -  V|),  (2.22) 

Cj  Cj 

where  S.  and  R’  indicate  the  positions  (z,y,z)  and  (2',y',2'),  the  horizontal  vector 
/  =  (i  —  x')i  +  (y  —  y')i  »  and  <7*  is  the  standard  deviation  of  the  random  quantity 
Scfci.  Further,  we  assume  that 

M(|z-z'|)dz'  =  2b,  (2.23) 

y-00 

where  Zq  is  a  measure  of  the  vertical  correlation  length,  which  is  assumed  much 
smaller  than  the  horizontal  correlation  length  Iq  given  below.  This  assumption 
is  based  on  the  fact  that  in  the  sediment,  the  medium  is  much  better  correlated 
horizontally  than  vertically.  Finally,  we  choose  the  correlation  fimction,  iV(I),  and 
its  associated  spectrum,  5(^),  to  have  one  of  the  following  forms: 


NiJ)  =  exp(-(|ri//o)], 

“  {1  +  (k/o)2F*’ 


{isotropic  exponential) 


(2.24) 


N{T) 

S(k) 


Ml 

lyO 


(1  +  {KxI^)^){1  +  («v^vo)*)’ 


{anisotropic  expor.ential) 


(2.25) 


n(T)  =  exp(— (/*//*o)’  —  (^v/^vo)*l>  {anisotropic  Gaussian) 

5(/e)  =  Jr/^/roexp[-^(«’/2o  +  «J^Jo))-  (2-26) 

Where  NiJ)  and  S{k)  are  related  as  Fourier  transforms: 
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S(k)  =  j  J  (2.27) 

For  the  isotropic  exponential  case,  Iq  is  the  correlation  length;  whereas  for  the 
anisotropic  exponential  case,  Igo  and  lyo  are  correlation  lengths  along  x  and  y  di¬ 
rections.  Even  when  l^o  =  lyo,  it  will  not  degenerate  to  the  isotropic  case,  since 
there  are  discontinuities  across  the  and  ly  axes.  The  discontinuity  problem  can 
be  corrected  by  introducing  a  second  length  8cale[l0],  however  this  will  not  be  dis¬ 
cussed  here.  As  for  the  Gaussian  case,  /xo.and  lyo  ue  the  same  as  in  the  anisotropic 
exponential  case.  However  when  /so  =  ^vOt  it  does  degenerate  to  the  isotropic  case. 
Since  anisotropic  correlations  are  Incorporated  in  the  formulation,  we  can  study 
the  anisotropic  effects  on  the  scattered  field,  which  will  be  covered  in  Chapter  4. 

We  notice  that  since  there  is  little  experimental  data  to  support  any  particular 
type  of  correlation  functions,  the  ones  we  choose  are  based  on  mathematical  sim¬ 
plicity.  The  correlation  length  in  all  three  types  of  correlation  fimctions  we  choose 
is  described  by  parameters  /o,  or  +  Z^.  Whether  the  modeling  is  enough  for 
real  ocean  bottoms  has  to  be  determined  by  geological  measurements.  However, 
the  formulation  developed  in  this  thesis  is  flexible  in  terms  of  fitting  other  types 
of  correlation  functions  as  long  as  they  are  given  in  the  form  of  two  dimensional 
spectrums. 

Now  we  have  completed  the  formulation  of  the  problem.  In  Chapter  3,  we 
will  solve  the  problem  of  coherent  reflection  when  a  plame  wave  is  incident  onto  the 
random  bottom,  and  in  Chapter  4  ,  we  will  concentrate  on  the  scattered  field.  Note 
that  we  ignored  the  second  order  term  in  £q.(2.10),  assuming  it  is  much  smaller 
than  the  first  order  term.  This  assumes  that  multiple  scattering  is  not  important. 
However,  when  the  sound  speed  variation  in  the  random  layer  is  so  large  that 
multiple  scattering  cannot  be  ignored,  the  present  approximation  will  fail. 
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Chapter  3 


COHERENT  REFLECTION 


In  Chapter  2,  we  decomposed  the  total  acoustic  wave  field  into  two  pairts:  the 
coherent  field  (P)  and  the  randomly  scattered  field  p*.  In  this  chapter,  we  will  use 
the  formulation  developed  in  the  previous  chapter  to  study  the  coherent  part  of 
the  total  field  (P)  when  a  plane  wave  in  an  oblique  angle  impinges  from  the  water 
to  the  bottom.  In  the  absence  of  the  random  layer,  the  refiection  problem  of  plane 
waves  from  the  interface  of  two  homogeneous  half-spaces  is  well  known.  As  intro¬ 
duced  in  Chapter  2,  we  have  a  random  layer  in  the  bottom  in  which  sound  speed 
has  a  random  component  (see  Fig.2.1).  Since  we  have  assumed  that  the  so\md 
speed  fluctuations  in  the  random  layer  are  small  compared  with  the  background 
bottom  sotmd  speed,  we  would  expect  that  there  is  a  small  loss  in  the  coherent 
component  of  the  sound  field. 

Assume  that  a  plane  wave  with  horizontal  wave  vector  kq  is  incident  upon 
the  bottom.  Since  the  coherent  pressure  field  in  the  water  obeys  the  Helmholtz 
equation,  it  has  to  take  on  the  follo^ndng  form: 

(Pi)  «  +  (K)e‘V*r^*)  ;  (3.1) 
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likewise,  the  coherent  wave  in  the  homogeneous  bottom  (medium  3)  will  have  the 
following  form: 

<Ps>  =  (r>e*‘®%-V^',  (3.2) 

where  {V)  and  (T)  are  the  coherent  reflection  and  coherent  transmission  coeffi¬ 
cients,  respectively.  They  are  constants  to  be  determined. 

Next  we  will  find  the  solution  of  (Ps) .  in  the  random  layer  and  then  use  the 
boundary  conditions  at  2  =  0  and  z  =■  —k  to  obtain  the  coherent  reflection  coeffi¬ 
cient  (V). 

Inserting  the  correlation  function,  Eq.(  2.22),  into  Eq.(  2.12),  the  integro- 
differential  equation  becomes: 

V>{P,{R))  +  k\{P,(R))  = 

~  j  1 1  -  z'm{S'))a,{R  \R')dv'.  (3.3) 

Since  the  incident  wave  is  in  the  kq  direction,  and  the  randomness  is  assumed 
stationary,  the  horizontal  wave  vector  Kq  in  the  coherent  wave  remains  constant 


throughout  the  medium.  Letting 

{P,{R))  = 

(3.4) 

we  obt^  from  Eq.(  3.3): 

(3.5) 

LI  \R')dx'd,/d2f. 


To  solve  for  P{z),  we  consider  the  case  Zq  <  A,  i.e.  the  acoustic  wavelength  is  much 
larger  than  the  vertical  correlation  length,  where  A  is  the  acoustic  wavelength. 
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Therefore,  we  can  write  the  vertical  correlation  aa  a  delta  function,  M{\z  —  z*!)  = 
zo^(2— 2'),  where  zq  is  a-  measure  of  the  vertical  correlation  length.  By  doing  so,  the 
integration  over  z'  is  eliminated.  The  assumption  that  the  random  sotmd  speed 
is  vertically  uncorrelated  is  supported  by  sediment  core  data[50].  Switching  the 
integration  variables  from  (x*,  y')  to  /  =  r  —  r*,  the  above  equation  can  be  simplified 
to: 

^  +  ti  -  (cjj  P(l)  = 

1 1  (3.6) 

Substituting  the  expression  for  G2  in  Eq.(  2.17)  into  Eq.(  3.6)  with  z*  -  z,  and 
rearranging  the  terms,  we  have: 

P(z)  =  -/(z)P(z)  (3.7) 


where  and  /(z)  are  defined  as  follows: 


T7’  =  kl(l-^+c(Ko)), 

(3.8) 

m  =  1  /  S(«- 

(3.9) 

Here  c(/Co)  is  a  constant  given  by 

(3.10) 

and  S{k  —  Kq)  is  the  shifted  two-dimensional  spectrum  of  N{l): 

S{iC-Kii)  =  1 1 

(3.U) 

Note  that  if  there  is  no  interface  discontinuity  at  z  =  0,  V21  in  Eq.(  3.9)  will  vauiish, 
likewise  /(z)  also  vanishes.  Eq.(  3.7)  then  becomes: 


This  is  a  homogeneous  Helmholtz  equation  with  complex  wavenumber  r).  There¬ 
fore,  in  the  absence  of  the  interface  discontinuity,  the  solutions  in  the  random  layer 
are  still  plane  waves,  but  with  complex  wavenumbers;  stated  in  another  way,  the 
solutions  a^-e  inhomogeneous  plane  waves.  Explained  physically,  the  random  layer 
subtracts  energy  from  the  coherent  propagation.  This  case  has  been  studied  by 
several  authors  mentioned  in  Chapter  1. 

By  including  the  water/bottom  interface  discontinuity,  we  have  an  extra  term, 
/(z),  in  Eq.(  3.7).  Now  plane  waves  are  no  longer  solutions  to  Eq.(  3.7)  since 
Eq.(  3.7)  is  no  longer  a  Helmholtz  equation  with  constant  wavenumbers.  However, 
observing  that  /(z)  is  always  much  smaller  than  since  Vn  is  a  quantity  whose 
magnitude  is  equal  or  smaller  than  one,  we  solve  Eq.(  3.7)  by  iteration.  The  zeroth 
order  solution  is 

Po(z)  =  yie-*"'  +  Be’”*,  (3.13) 

and  the  first  order  solution  is 

P(z)  =  +  Be’"*  -  ^  f{z’)Po{z')e^’’\‘-’’\dz' 

=  +  Be’"*  -  ^  /(z')(yle-’'"'  +  Bc^’^‘')e'^^—''>dz’ 

r  /(z')(Ae-"»''  -h  Be’^^-’^-'Jdz',  (3.14) 

2tT]  J 

where  A  and  B  are  unknown  constants.  The  integration  over  depth  z'  will  be  car¬ 
ried  out  numerically. 

Then  the  approximate  solution  of  the  coherent  pressme  field  in  the  random 
layer  is: 

{P,{R))  =  P{2)e<'^'.  (3.15) 
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Now  we  have  solutions  for  the  coherent  fields  in  media  1,  2  and  3,  i.e.,  Eqs.(  3.1), 
(  3.2)  and  (  3.15),  with  four  unknown  constants  (V),  (T),  A  and  B.  By  requiring 
the  continuity  of  the  mean  pressure  and  mean  vertical  particle  velocity  along  the 
interfaces  z  =  0  and  z  a=  — h. 


(P.(r,0))  =  (i»,(nO)), 

(3.16) 

(3.17) 

(P,(r,-k))  =  <Ps(r,-A)), 

(3.18) 

(3.19) 

we  obtain  the  coherent  refiection  coefficient  (V): 

(*>01  +  Vll  +  r6Voit>23«**'’*  +  <12V2S«’*'''‘) 

^  (1  +  VolUlJ  +  +  txjt?0lV23C**'''*)  ’ 

(3.20) 

In  the  above  expression,  there  are  several  quantities  that  will  be  defined.  Most  of 

them  involve  the  quantities  q  or  <r  due  to  the  small  fiuctuations  of  soimd  speed  in 

the  layer.  They  are: 
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j  (\/fc|-ic§  +  »7)a 

(3.25) 

The  quantities  a,  7,  (  and  &re: 

(3.26) 

(3.27) 

7  =  1-  f{z')e^^^'^dz\ 

(3.28) 

(3.29) 

2iri  J-h 

(3.30) 

Recall  that  the  quantity  f{z)  is  given  by  Eq.(  3.9). 

By  inspection,  we  observe  that  when  a  approaches  zero,  c(k^)  and  f[z)  also 
approach  zero,  17  =  yjk\  —  Kqi  wid  ^  become  unity,  {P2[R))  degen¬ 

erates  to  the  conventional  transmitted  wave,  and  (V)  is  merely  the  conventional 
plane-wave  reflection  coefficient  for  two  half-spaces. 

The  coherent  reflection  coefficient,  (V),  in  Eq.(  3.20)  is  obtained  by  using 
Eqs.(  3.21)  through  (3.30)  and  Eq.(  3.9)  through  numerical  integrations.  The  mod¬ 
ulus  and  phase  of  (V)  versus  incident  angle  9q  for  difierent  values  of  the  statistical 
parameters,  including  the  case  of  non-random  reflection  for  comparison,  are  shown 
in  Figs.3.1-3.6.  Among  them,  Figs.3.1-3.4  are  for  the  case  of  a  Tast”  bottom,  i.e. 
C]  >  Cl,  and  Fig8.3.5-3.6  are  for  the  case  of  a  “slow**  bottom,  Cj  <  ci.  Figs.3.1 
and  3.2  are  the  magnitude  and  phase  of  the  coherent  reflection  coefficient  for  three 
different  frequencies,  200  Hz,  300  Hz  and  400  Hz.  The  top  line  is  the  reflection 
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coefficient  with  no  randomness.  The  results  show  an  increased  energy  loss  with 
increasing  frequency.  Figs.3.3  and  3.4  show  the  magnitude  and  phase  of  the  coher¬ 
ent  reflection  coefficient  for  the  same  set  of  acoustic  background  parameters  with 
flxed  frequency,  but  changing  randomness  size,  respectively,  o  =  0,<t  s  0.05  and 
a  —  0.10.  Here  we  see  the  energy  loss  increases  with  increasing  a.  Compared  with 
the  case  of  reflection  from  a  fast  homogeneous  half-space,  the  major  effect  due  to 
the  randomness  is  the  disappearance  of  a  sharp  critical  angle.  We  also  notice  that 
the  phase  change  is  not  significant.  Figs.3.5  zmd  3.6  are  the  magnitude  and  phase 
of  (V)  versus  different  frequencies  for  a  slow  bottom.  In  this  case,  total  reflection 
does  not  occur.  The  change  in  magnitude  is  not  as  significant  as  for  fast  bottoms. 
However,  the  phase  change  is  relatively  large. 

The  behavior  of  (V)  is  similar  to  the  reflection  coefficient  for  a  homogeneous 
absorbing  bottom[56].  However,  the  mechanism  of  energy  loss  is  different[58]. 

After  finding  (V),  the  other  unknowns.  A,  B,  and  (T),  can  be  found  through 
Eqs.(  3.16)  to  (  3.19)  to  be: 

^  -  P2v/fci  -«o));  (3.31) 

B  =  (3.32) 

(2^  =  expl-i^k^  -  K^h]laAexp(trjh)  +  7Bexp(-irjh)].  (3.33) 

Notice  that  when  <7  approaches  zero,  (T)  degenerates  to  the  conventional  transmis¬ 
sion  coefficient.  These  results  will  be  ;ised  in  Chapter  4  to  calculate  the  scattered 
field. 
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{a)  iO  aSVHd 
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Figure  3.2:  Phase  of  reRection  coeflicient  for  difTerent  frequencies  (fast  bottom) 
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Figure  3.3:  Magnitude  of  reflection  coelRcient  for  different  a'a  (fast  bottom) 


(a)  jo  asvHd 
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Figure  3.4:  Phase  of  reflection  coeflicient  for  different  o’s  (fast  bottom) 


42 


Figure  3.5:  Magnitude  of  reflection  coeflicient  for  different  frequencies  (slow  bot- 


{a)  dO  aSVHd 
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Figure  3.6:  Phase  of  reflection  coeflicieni  for  different  frequencies  (slow  hot- 


Chapter  4 


SPATIAL  CORRELATION  OF 
THE  SCATTERED  FIELD 


In  Chapters  2  and  3,  we  have  found  the  solutions  of  the  coherent  reflection 
and  transmission  coefficients  {V)  and  (T),  and  the  constants  A  and  B,  hence  the 
coherent  fleld  {P2{M'))  is  known.  We  can  put  the  coherent  fleld  into  £q.(  2.11)  in 
Chapter  2  to  find  the  randomly  scattered  field: 

p:(^) ix')dv'.  (4.i) 

We  are  interested  in  the  scattered  acoustic  field  in  the  water,  as  measured  at  the 
observation  point  R  which  is  in  the  water  column  where  z  >  0.  The  only  difference 
between  Eq.(  2.11)  and  Eq.(  4.1)  is  that  the  observation  point  in  Eq.(  2.11)  is  in 
the  random  layer  where  z  <  0.  Therefore,  the  Green’s  function  that  will  be  used 
here  is  one  with  the  source  in  the  random  layer  {z  <  0)  and  the  observation  point 
in  the  water  (z>0). 

Since  the  statistical  average,  i.e.  the  first  moment,  of  the  random  quantity  e  is 
zero,  (e)  =  0,  (p*)  =  0  becaxise  p*  is  proportional  to  c.  Hence  we  shall  study  the 
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sp&tial  correlation  function,  i.e.  the  second  moment  of  p*,  which  is  defined  as: 

caSuS,)  =  o>;(^.)pr(^)>.  (■‘■2) 

where  Ri  and  Rt  are  two  spatial  positions  where  the  receivers  are  located.  There¬ 
fore,  C«(^i,^2)  is  &  measure  of  how  well  the  scattered  fields  at  different  loca¬ 
tions  are  correlated.  Notice  that  the  the  spatial  correlation  of  the  scattered  field, 
C«(Ai,J?2),  defined  above  is  not  normalized;  when  Ri  =  R3,  is  the  scattering 
intensity. 

Conventionally,  only  C^{Ri,Ri)  —  (|p*(i2i){*),  the  intensity  of  the  scattered 
field,  is  studied.  However,  the  intensity  contains  just  a  small  portion  of  the  in¬ 
formation  carried  in  the  second  moment  of  the  scattered  field.  No  effort  has  yet 
been  made  to  study  the  spatial  correlation  of  the  quantity  p*  in  bottom  acoustics. 
Results  of  this  kind  exist  in  electromagnetic  wave  and  atmospheric  acotistic  wave 
studies[16],  where  the  media  are  assumed  infinite,  and  no  medium  discontinuities 
are  present.  The  formulation  given  in  Chapter  2  allows  us  to  calculate  the  spatial 
correlation  of  the  scattered  field,  not  only  the  intensity.  Recent  technological  de¬ 
velopments  in  underwater  acoustic  experiments  make  it  possible  to  investigate  the 
spatial  correlation  of  the  scattered  field.  Some  considerations  on  experiments  of 
this  kind  will  be  given  in  Chapter  7. 

Substituting  Eq.(  4.1)  into  Eq.(  4.2),  we  have 

caSuS,)  =  [%)'!  a 

{P3{R')){p;{R”))g{Rx\r)g-{R2\r”).  (4.3) 

The  correlation  function  of  e  is  given  in  Eq.(  2.22)  of  Chapter  2.  If  the  acoustic 
wavelength  is  much  larger  than  the  vertical  correlation  length  of  e,  A  >>  sq*  we 
can  approximate  the  correlation  function  as  follows: 
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(€(5)£(fi')>  =  ^  AzQa^N{T)6{z  -  z'),  (4.4) 

C2  Cj 

where  N(t)  and  its  spectrum  are  assumed  to  have  one  of  the  forms  given  in 
Eqs.(  2.24)  to  (  2.26).  The  Green’s  fimction  in  which  the  source  is  located  in 
the  random  layer  and  the  observation  point  is  located  in  the  water  has  this  form: 

i  /  /»(«,*!,  (4.5) 

where  p(/?,  Zi,  s')  is  the  two  dimensional  Fourier  transform  of  the  Green’s  Fimction, 
G(iZi|i2'),  and  is  sometimes  called  the  depth-dependent  Green’s  function[54].  For 
the  present  case,  it  is: 

,(«,  »-)  =  r„  f  (4.6) 

where  (x  uid  ^3  are  vertical  wavenumbers  in  the  water  and  in  the  bottom,  respec¬ 
tively  corresponding  to  k  : 

fa  =  y/kl-  K*,  (4.7) 


and  T21  is  the  transmission  coefficient  from  medium  2  to  medium  1  in  the  absence 
of  the  the  random  quantity  e: 


r,i  = 


2pifj 


Pifa  +  Pafi 

Substituting  the  above  expressions  into  Eq.(  4.3),  we  have 


(4.8) 


5  («,  *1,  Z2,  s')e’*  (4.9) 
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The  above  formula  is  the  basis  for  the  calculation  of  the  spatial  correlation  of  the 
scattered  field.  It  consists  of  two  wavenumber  integrals  and  one  depth  integral 
over  the  remdom  layer.  1P(^')I  z-dependent  part  of  the  coherent  field  in  the 
random  layer;  it  can  be  calc\ilated  through  Eqs.(  3.14),  (  3.31),  (  3.32),  and  (  3.33) 
in  Chapter  3.  As  a  comparison,  we  can  also  use  the  Extended  Bom  approximation 
in  Eq.(  4.9).  To  do  so,  instead  of  using  the  coherent  solution,  we  use  the  backgroimd 
field,  i.e.  the  field  in  the  absence  of  the  randomness: 

P{z)  =  rx3e-*V*^*,  (4.10) 


where  Tu  is  the  transmission  coefficient  from  medium  1  to  2: 

Later  in  this  chapter,  a  comparison  between  using  the  coherent  field  and  using  the 
Extended  Bora  approximation  will  be  made. 


In  order  to  integrate  Eq.(  4.9)  numerically,  we  chamge  the  wavenumber  integrals 
into  polar  coordinates.  Letting  the  polar  coordinates  for  the  quantities  ic,  acq  and 
r  —  respectively  be: 

K=  (k,^), 

Ko  =  («o,^o), 

f-ro  =  {R,e),  (4.12) 

Eq.(  4.9)  becomes: 

J  d<pS{Kcos<^  -  Kocos^^KStTKp  -  (4.13) 

where  R  is  the  horizontal  separation  distance  of  the  two  receivers,  and  0  is  the  hori¬ 
zontal  angle  between  the  line  linking  the  two  receivers  and  the  x-axis.  The  following 
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results  for  the  spatial  correlation  of  the  scattered  field  are  obtained  from  Eq.(  4.13). 

Fig.4.1  shows  the  receiver  geometry.  Since  the  two  receiver  points  in  Eq.(  4.9) 
are  completely  arbitrairy,  we  will  choose  them  to  be  either  in  a  vertical  line  parallel 
to  the  z-axis,  or  a  horizontal  line,  perpendicular  to  the  z-axis.  The  orientation 
of  the  horizontal  receivers  is  indicated  by  the  azimuthad  angle  B.  The  incident 
wave  direction  is  determined  by  the  incident  angle  6q  and  the  azimuthal  angle 
^0*  Although  the  anisotropy  of  the  bottom  correlation  function  is  assumed  to  be 
aligned  in  the  xy-plane  (see  Eqs.(  2.24)  to  ( 2.26)),  neither  the  incident  wave  nor  the 
horizontal  receivers  have  to  be  so  aligned;  therefore  this  configuration  is  arbitrary 
in  terms  of  the  incident  wave  direction  and  receiver  orientation. 

For  clarity,  Fig.4.2  specifies  the  environmental  parzLmeters  upon  which  all  the 
following  results  are  calculated  except  for  those  separately  indicated.  Since  there  is 
little  experimental  data  available,  especially  those  concerning  the  random  distribu¬ 
tion  of  sound  speed  in  the  layer,  the  parameters  are  chosen  to  meet  the  theoretical 
requirements  and  to  be  in  the  reasonable  range  of  reality.  How  to  experimentally 
measure  these  parameters  will  be  discussed  in  Chapter  7. 
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Fig.  4.1  Incident  wnve  and  receiver  axrangexnent 
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Fig.  4.2  Environment  parameters 


4.1  Angular  Dependence 


First,  we  examine  the  angular  dependence  of  the  horizontal  correlation.  In  the 
present  case,  we  choose  the  incident  angle  9o  to  be  either  30**  or  60^,  and  =  0,  so 
the  incident  wave  is  in  the  x  direction.  The  bottom  correlation  function  is  isotropic 
Gaussian  (see  Eq.(  2.26)).  In  calculating  the  horizontal  correlation  of  the  scattered 
field,  fi),  we  keep  one  receiver  at  a  fixed  point,  and  the  second  receiver  takes 
sequential  meastirements  away  from  the  first  receiver  in  the  horizontal  plane.  The 
distance  between  the  two  receivers  is  R.  The  second  receiver  may  also  move  in 
different  directions,  which  correspond  to  different  orientation  angles.  Since  the 
incident  direction  is  fixed  in  the  x-direction,  we  are  going  to  examine  the  differ¬ 
ences  in  horizontal  correlations  among  different  receiver  orientations.  Throughout 
Fig.4.3  to  Fig.4.8,  the  acoustic  frequency  is  100  Hz.  The  x-direction  is  referred 
to  as  along-propagation  direction,  and  the  y-direction  is  referred  to  as  the  cross- 
propagation  direction.  Fig.4.3  shows  the  comparison  between  the  horizontal  cor¬ 
relation  in  the  along-propagation  direction  (0  =  0®)  and  the  horizontal  correlation 
in  the  cross-propagation  direction  (0  =  90®)  while  the  incident  angle  =  30®.  The 
horizontal  axis  is  expressed  in  the  non-dimensional  unit  R/X.  Notice  that  while 
the  cross-propagation  correlation  (0  =  90®)  has  a  clear  null  near  0.6  wavelengths, 
the  along-propagation  correlation  (0  =  0®)  does  not.  In  other  words,  the  corre¬ 
lation  along  the  propagation  direction  is  greater  than  the  correlation  across  the 
propagation  direction.  For  better  views,  Figs.4.4  and  4.5  show  two  side  views  of 
the  gradual  change  of  |Cv(i2,9)|  against  different  9's  ranging  from  G®  to  90®.  One 
of  the  two  axes  is  still  R/X,  the  other  one  is  0. 

Figs.4.6,  4.7  and  4.8  are  the  same  as  Figs.4.3,  4.4  and  4.5  except  that  the  inci¬ 
dent  amgle  is  increased  to  60®. 
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Cheniov[l6]  mentioned  a  similar  angular  dependence  for  the  case  of  scattering 
in  infinite  media.  Physically,  since  the  incident  wave  is  a  plane  wave,  receivers 
across  the  propagation  direction  will  have  the  same  phase,  whereas  along  other 
directions  the  phase  differs,  and  the  difference  reaches  maximum  when  the  receivers 
are  aligned  along  the  propagation  direction.  Later  it  will  be  shown  that  the  null  in 
the  correlation  can  be  used  to  estimate  the  bottom  correlation  length  by  an  inverse 
procedure.  Therefore  in  an  experiment,  the  receivers  should  be  deployed  across  the 
incident  wave  propagation  direction  to  identify  the  null. 
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Comparison  on  receiver  mientation 


Magnitude  of  C(R,  6  ) 


Figure  4.4:  Comparuon  on  receiver  orientation,  fio  s  30* 
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f  «  100  Hz 


Figure  4.5:  Comparison  on  receiver  orientation,  tfo  *  30* 
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Comparison  on  receiver  orientation,  6o  — 


Magnitude  of  C(R,  6  ) 


Figure  4.8:  Comparison 


on  receiver  orientation,  =  60* 


4.2  Frequency  Dependence 


Next  we  investigate  the  frequency  dependence  of  the  spatial  correlation  func¬ 
tion  Cv{Ri9).  We  will  emphasize  the  relationship  between  the  physical  size  of 
the  scatterer  expressed  by  twice  the  horizontal  correlation  length,  i.e.  2/o,  and  the 
acoustic  wavelength  used.  Again  we  choose  the  correlation  function  to  be  Gaussian 
in  this  section.  Since  we  have  assumed  that  =  S  meters  throughout  this  section, 
the  physical  size  of  the  volume  scatterer  has  a  diameter  of  roughly  10  meters.  For 
comparison,  we  vary  the  acoustic  frequencies.  Because  the  acoustic  wavelength  A  is 
a  function  of  frequency,  it  is  not  convenient  to  use  the  unit  RJX  for  the  horizontal 
axis  as  we  did  in  the  last  section;  therefore,  the  horizontal  axes  of  the  Figures  in 
this  section  are  all  expressed  in  terms  of  iZ  in  the  units  of  meters. 

Fig.4.9  shows  the  horizontal  correlations  when  the  incident  angle  9q  is  30®  for 
three  different  frequencies,  40  Hz,  60  Hz  and  80Hz,  respectively.  The  correspond¬ 
ing  acoustic  wavelengths  are  37.5m,  25m  and  18.75m  respectively  and  they  are  all 
larger  than  the  volume  scatterer  size  of  10  meters.  We  find  that  the  first  null  of  the 
correlation  changes  accordingly  when  the  acoustic  frequency  increases.  It  changes 
in  such  a  fashion  that  the  null  position  on  the  RfX  axis  will  be  fixed.  In  Fig.4.10, 
the  acoustic  frequencies  are  160  Hz,  180  Hz  and  200  Hz,  and  the  corresponding 
wavelengths  are  9.38m,  8.33m  and  7.5m  respectively.  They  are  all  smaller  than 
the  scatterer  size  of  10  meters.  This  time  the  null  positions  are  independent  of 
frequency;  more  importantly,  the  half-value  point  on  the  correlation  curve  roughly 
equals  to  5  meters,  which  is  the  correlation  length  /q  of  the  random  sound  speed. 

These  results  show  that  when  the  acoustic  wavelength  is  comparable  or  smaller 
than  the  size  of  the  scatterer,  the  position  of  the  first  null  is  independent  of  fre- 
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quency,  and  the  half-value  point  equals  the  correlation  length  of  that  of  e.  However, 
when  the  acoxxstic  wavelength  is  larger  than  the  size  of  the  scatterer,  the  first  null 
changes  with  frequency. 

In  order  to  see  this  more  clearly,  we  show  in  Fig.4.11  a  three  dimensional  plot 
of  the  correlation  Cv{R,6).  The.  two  horizontal  axes  are  correlation  distance  R  in 
meters  and  frequency  in  Hz.  Note  the  two  regions  where  the  first  null  locus  has 
different  behavior  against  frequency;  it  curves  into  smaller  values  for  increasing  fre¬ 
quency  starting  from  the  low  frequency  end;  and  as  the  frequency  increases  further 
when  the  acoustic  wavelength  is  smaller  than  10  meters,  the  null  position  stablizes 
at  10  meters. 

Physically,  when  the  acoustic  wavelength  is  larger  than  the  linear  dimension  of 
the  scatterer,  the  wave  cannot  distinguish  the  scatterer  size,  and  so  what  is  reflected 
on  the  correlation  of  the  scattered  field  of  the  random  scatterers  is  the  character 
equivalent  to  what  of  a  white  random  process[65]  would  have.  When  the  acoustic 
wavelength  becomes  smaller  than  the  linear  dimension  of  the  volume  scatterer,  the 
spatial  correlation  of  the  scattered  field  will  reflect  the  size  of  the  volume  scatter. 

Figs.4.12,  4.13,  and  4.14  are  the  same  as  Figs.4.9,  4.10,  and  4.11  except  that 
the  incident  angle  is  60”  rather  than  30”.  Figs.4.15,  4.16,  and  4.17  are  also  the 
same  as  Figs.4.9,  4.10,  and  4.11,  but  show  the  vertical  correlation,  rather  than 
the  horizontal  correlation.  The  label  Z  in  the  figures  is  the  vertical  separation 
distance  between  the  two  receivers  in  meters.  In  all  these  figures,  the  characteristic 
frequency  dependence  of  the  correlation  is  the  same.  This  is  a  very  important 
result  in  terms  of  inverting  the  ocean  bottom  parameters.  If  we  sweep  the  acoustic 
frequency  in  an  experiment  to  measure  the  spatial  correlation  of  the  scattered  field, 
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and  form  the  correlations  on  either  horizontal  arrays  or  vertical  arrays  for  different 
frequency  bins,  >ve  should  find  a  dividing  point  at  which  the  first  null  will  stop 
changing  with  increasing  frequency.  That  point  will  indicate  the  correlation  length 
of  the  bottom  random  sound  speed  distribution.  In  a  real  experiment,  there  will 
be  ambient  noise  compoimded  into  the  scattering  field  and  make  the  identification 
of  the  nulls  difficult.  However,  we  can  always  pick  the  half-value  point  on  the 
correlation  curve  to  estimate  the  size  of  the  scatterers. 
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Figure  4.11:  Frequency  dependence,  —  30® 
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Figure  4.12:  Frequency  dependence,  low  frequency  case,  tfo  =  ®0* 
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Frequency  dependence,  high  frequency  case,  io  ~ 
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Frequency  dependence,  low  frequency  case, 
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Figure  4.16;  Frequency  dependence,  high  frequency  case,  =  30* 


Figure  4.17;  Frequency  dependence,  $q  =  30* 


4.3  Comparison  with  the  Born  approximation 


The  Bom  approximation  has  been  used  extensively  in  the  calculation  of  weak 
scattering  problems.  However,  when  the  size  of  the  scatterer  is  large  compared  to 
the  wavelength,  the  Bom  approximation  fails  since  it  does  not  include  the  energy 
loss  in  the  incident  wave  due  to  scattering.  The  bi-local  approximation  used  in 
this  thesis  takes  into  account  the  loss  of  energy  in  the  coherent  part  of  the  field  as 
discussed  in  Chapters  2  and  3.  In  this  section  we  compzire  the  two  methods.  In 
£q.(  4.13),  if  the  bi-local  approximation  is  tised,  the  quantity  |P(z')|  is  given  by 
Eq.(  3.14]  in  Chapter  3;  if  the  Extended  Born  approximation  is  used,  |P(s')|  is  given 
by  Eq.(  4.10]  in  this  chapter.  The  difference  between  the  two  approximations  is 
that  the  vertical  wavenumber  in  the  Bom  approximation  is  real,  and  therefore  there 
is  no  energy  loss  while  propagating  deep  into  the  random  layer,  whereas  the  vertical 
wavenumber  is  complex  in  the  bi-local  approximation,  hence  the  wave  is  evanescent, 
and  wave  energy  is  taken  away  due  to  scattering  while  propagating  deep  into  the 
random  layer.  Fig.4.18  shows  the  horizontal  correlations  of  the  scattered  field  when 
the  thickness  of  the  random  layer  equals  six  acoustic  wavelengths.  We  see  little 
difference  between  the  two  results.  In  Fig.4.19,  all  parameters  are  kept  the  same  as 
in  Fig.4.18  except  that  the  layer  thickness  b  sixty  times  the  acoustic  wavelength. 
There  we  find  the  horizontal  correlation  of  the  scattered  field  b  apparently  different. 
The  Born  approximation  overestimates  the  scattering  intensity  (RsO).  At  the  same 
time  we  notice  that  the  correlation  patterns  obtained  from  the  two  methods  are 
essentially  the  same.  Figs.4.20  and  4.21  repeat  the  results  in  Figs.4.18  and  4.19, 
except  they  are  vertical  correlations,  and  likewbe  we  find  the  results  are  similar. 
Therefore,  we  can  conclude  that  the  Bom  approximation  can  be  used  as  a  fast 
way  to  examine  the  correlation  pattern,  whereas  when  the  scattering  intensity  b 
concerned,  the  Bom  approximation  will  result  in  more  error.  If  the  random  layer 
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is  a  random  half-space,  the  Bom  approximation  will  not  converge,  but  the  bi-local 
approximation  will  still  be  applicable. 
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Figure  4.18:  Comparison  with  the  Bom  approximation,  thin  layer,  horisontal  cor¬ 
relation 
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Figure  4.19:  Comparison  with  the  Born  approximation,  thick  layer,  horizontal 
correlation 
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Figure  4.20:  Comparison  with  the  Born  approximation,  thin  layer,  vertical  corre¬ 
lation 
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Comparison  with  the  Born  approximation,  thick  layer,  vertical  corre> 


4.4  The  Anisotropy 


So  far  we  have  only  studied  the  case  where  the  bottom  correlation  functions  is 
isotropic,  i.e.  lto  =  lyo-  It  is  also  important  to  examine  the  influence  of  anisotropic 
correlation  of  e  on  the  spatial  correlation  of  the  scattered  field.  The  formulsi- 
tion  developed  in  this  thesis  allows  such  anisotropic  correlations.  In  the  following 
examples,  the  environmental  parameters  will  be  kept  the  same  as  in  Fig.4.2  except: 

1)  Igo  =  2.0m,  and  lyo  =  10.0m,  so  the  sound  speed  variation  in  the  random 
layer  is  more  correlated  in  the  y-direction  than  in  the  x-direction. 

2)  In  Figs.4.22-4.24,  the  correlation  functions  are  still  Gaussian,  but  in  Figs.4.25- 
4.27,  the  correlation  functions  are  2-d  exponential  type. 

All  the  examples  are  horizontal  correlations  of  the  scattered  field  measured 
across  the  propagation  direction. 

Figs.4.22,  4.23,  and  4.24  show  the  examples  for  incident  wave  frequencies  of 
50,  100,  and  200  Hz,  respectively.  The  curves  corresponding  to  ^  =  0  are  those 
when  the  incident  wave  is  in  the  x-direction,  therefore  confronting  the  wide  side  of 
the  volume  scatterer;  the  other  curves  corresponding  to  =  90*’  are  those  when 
the  incident  wave  is  in  the  y-direction,  confronting  the  narrow  side  of  the  volume 
scatterer.  We  see  that  the  ^  =  0  curves  all  indicate  a  10  meter  correlation  length 
on  their  first  nulls.  It  is  clear  that  when  the  incident  wave  comes  in  the  direction 
of  the  wide  side  of  the  scatterer,  the  longer  correlation  length  will  be  reflected  on 
C„{R,9).  We  notice  that  the  nulls  on  the  ^  =  90**  curves  change  with  frequency. 
This  is  because  when  the  incident  wave  is  in  the  direction  of  the  narrow  side  of  the 
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scatterer,  the  null  position  will  not  stablize  until  the  acoustic  wavelength  is  shorter 
than  the  smaller  correlation  length,  i.e.  2  meters  in  these  examples.  Figs.4.25, 
4.26,  and  4.27  are  the  same  as  Fig.4.22,  4.23,  and  4.24,  except  that  the  Gaussian 
correlation  function  is  changed  to  2-d  exponential  correlation  function.  We  find 
similar  results  in  these  examples. 

We  conclude  that  anisotropic  correlation  functions  will  infiuence  the  spatial 
correlation  functions  of  the  scattered  field  considerably.  When  the  null  position  is 
used  to  invert  for  the  correlation  length,  it  is  necessary  to  change  the  direction  of 
the  incident  wave  to  see  whether  anisotropy  exists. 
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Figure  4.22:  Anisotropy  effect,  f=50  H* 
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Figure  4.26:  Anisotropy  effect,  f=100  H* 
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Anisotropy  effect,  f=200  He 


In  Chapters  2,  3,  and  4,  we  studied  the  wave  scattering  from  a  random 
bottom.  One  major  goal  in  this  research  is  to  invert  for  the  bottom  sound  speed 
parameters.  We  concluded  that  the  correlation  length  of  the  bottom  random  soimd 
speed  is  related  to  the  spatial  correlation  length  of  the  scattered  field.  It  b  worth¬ 
while  to  compare  this  research  with  work  on  scattering  from  a  deterministic  inho¬ 
mogeneous  ocean  bottom.  There  b  an  enormoTis  literature  on  inverse  problems  of 
thb  kind.  Notably,  Cohen  and  Bleinstein[59,60,61,62],  along  with  other  researchers, 
studied  the  velocity  inversion  problem  when  the  velocity  variation  b  non-random 
and  small.  They  expressed  the  soxmd  velocity  variation  through  an  integral  equa¬ 
tion  related  to  received  signal  by  using  the  Bom  approximation.  The  studies  of 
random  scattering  problems  can  be  linked  to  the  studies  of  non-random  scattering 
problems  by  introducing  randonmess  to  the  formulated  non-random  problem  and 
properly  averaging  the  results  to  invert  for  the  bottom  parameters  related  to  the 
non-random  problem. 
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Chapter  5 


VOLUME  SCATTERING  VS. 
ROUGH  INTERFACE 
SCATTERING 


As  mentioned  in  the  Introduction,  both  bottom  inhomogeneities  and  interface 
roughness  contribute  to  ocean  bottom  scattering.  If  the  interface  is  very  rough, 
and  has  length  dimensions  which  are  large  relative  to  the  acoustic  wavelength, 
one  can  expect  that  rough  interface  scattering  will  dominate  the  scattered  field; 
if  the  roughness  is  small  relative  to  the  acoustic  wavelength  and  gentle  in  slope, 
the  scattered  field  may  be  the  combined  result  of  volume  and  interface  scattering 
mechanisms. 

Since  there  is  no  existing  theory  that  models  both  the  roughness  and  the  volume 
inhomogeneity  in  one  unified  approach,  there  is  little  quantitative  understanding  of 
the  relative  importance  of  the  two  scattering  mechanisms.  This  chapter  is  intended 
to  Investigate  the  relative  strengths  of  the  scattered  fields  resulting  from  interface 
roughness  and  from  volume  inhomogeneities,  and  if  both  scattering  mechamisms 
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Figure  5.1:  Geometry  of  roughness  on  top  of  volume  inhomogeneity 

are  important  in  one  location,  whether  it  is  possible  to  distinguish  between  the  two 
scattering  mechanisms. 

The  geometry  is  shown  in  Fig.5.1.  We  assume  that  the  water/bottom  interface 
roughness  is  random  and  small  in  amplitude  and  slope,  and  the  bottom  sound 
speed  variation  Sc  is  random  and  small  and  superimposed  upon  a  constant  as 
in  the  previous  chapters.  Further  we  assume  that  there  are  no  statistical  correla¬ 
tions  between  the  interface  roughness  and  the  volume  inhomogeneities.  The  total 
randomly  scattered  field  can  be  written  as  the  summation  of  the  interface 
scattered  field  p*  and  volume  scattered  field  p*: 

pUui  *  pI  +  PJ-  (5.1) 

Therefore,  the  second  moment  of  the  total  scattered  field  is: 

<pr«-ps..i)  -  <(p:+p:)(pr+pr)> 

-  w:-^j^.k+kpv+p':p‘.) 
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Figure  5.2:  Schematic  diagram  of  wave  scattering  from  interface 

=  w)  +  «?:•).  (5-2) 

where  the  cross  correlation  terms  between  the  interface  scattered  field  and  the 
volume  SC'  iered  field  have  been  omitted  because  the  two  random  processes  are 
assumed  statistically  independent.  So,  the  correlation  of  the  total  scattered  field 
is  simply  the  summation  of  the  correlations  of  the  interface  scattered  field  and  the 
volume  scattered  field.  This  means  that  the  two  scattering  mechanisms  can  be 
examined  separately  under  the  above  mentioned  assumptions. 

The  Bass  and  Fuks’[l9]  perturbation  theory  will  be  used  to  find  the  scattered 
field  correlation  caused  by  the  rough  interface  alone,  and  the  spatial  correlation  of 
the  volume  scattered  field  obtamed  in  the  previous  chapters  will  be  combined,  to 
investigate  the  combined  effects. 

According  to  Bass  and  Fuks[19|,  let  two  media  be  separated  by  the  rough  in¬ 
terface  s  =  y)  sound  speed  and  density  in  each  meditim  are  given  in 
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Fig.5.2.  If  we  designate  P,  as  the  acoustic  pressure,  the  vertical  particle  velocity 
and  the  pressure  are  continuous  across  the  interface,  t.e. 

where  n  is  the  normal  of  the  interface.  We  decompose  P,i  and  P,2  into  the  following 
forms: 

P.i  =  P!i  +  P.i;  Pa  =  Jf,  +  Pa,  (5-4) 

where  P^i  and  are  the  pressure  in  the  upper  and  lower  media  in  the  absence 
of  the  roughness,  while  p^i  and  p«2  the  components  related  to  the  presence  of 
small  roughness. 


Taking  into  account  that 

and  noticing  that  along  z  s=  0, 

7^~dr  -  TrTaT' 

we  expand  the  boundary  conditions  of  £q.(  5.3)  as  a  power  series  in  ^  and  7,  and 
retain  the  terms  up  to  first  order.  Then  the  boundary  conditions  found  for  the 
scattered  fields  p«i  and  p,i  at  z  =  0  are 

P.i  “  P.a  =  p{f) 

where  p{f)  and  v(f),  the  pressure  and  vertical  particle  velocity  discontinuities, 
respectively,  along  z  =  0,  are 


P1P2  Pi-'  Pi  , 
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and  ki  and  are  the  wavenumbers  in  the  upper  and  lower  media,  respectively, 

If  the  scattered  field  is  written  in  the  form  of  plane  waves  leaving  the  averaged 
interface,  z  3=  0: 

p,i{R)  =  [  =  \Jkl-K^ 

00 

P»7[R)  =  f  /Cj,  =  y/kl-K?  (5.9) 

y-00 

we  then  apply  the  Fourier  transform  to  both  sides  of  Eq.(  5.7): 


p;i(«)  -  Prf(«)  =  P(^), 

+  ^^p72(«)  =  5(^, 

Pi  Pj 

and  note  that  p«i(k)  can  be  expressed  as 

PiK2.p(«)  -  *PiP25(k) 

P*i('cJ  —  V  -  "  » 

PtKit  -h  PiKu 

where  p  and  v  are  Fourier  transforms  for  p{f)  and  v(f),  respectively. 


(5.10) 

(5.11) 


Let  a  plane  wave  in  the  (Jb^,  kyo,  -kgo)  direction  be  the  incident  wave,  then 
pO^  _  g«(*»o*+*fOV-*«o»)  ^  pjjg«(*.o*+ktO»+*.o»)^  (5*12) 


and  is  the  plane  wave  refiection  coefficient: 

V„  =  (S.13) 

Pik^  +  PiV^  -kl  +  klo 

Substituting  £q.(  5.12)  into  Eq.(  5.8),  and  applying  the  Fourier  transform  to  both 
sides,  we  find  that  p  and  v  take  on  the  following  forms: 

p(«)  =  -ik^{l-Vo)^^^^^n'^,-k^,Ky-kyo) 

Pi 

5(;c)  =  (i  +  Vb)i^2J^?(K.-k^,S-A:,o) 

P1P2 

-kgOKg  —  kyOKy  +  — !h£i  .  (5.14) 

P2-'Pl  J 

In  these  expressions,  ^  is  the  Fourier  transform  of 


(5.14) 
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(5.15) 


Now,  p',i  in  Eq.(  5.11)  can  be  solved  by  using  Eq.(  5.14) 

Pii  =  ~~i — ~  kzOti^v  ~ 

PiKu  +  Pl'Cj* 

•  -  Vq)  +  (1  +  Va)i-kxoKx  -  kyoKy  +  _  p*^"]  * 

The  scattered  field  in  the  upper  medium  is  the  inverse  Fourier  transform  of  p,i: 

Pti(^)  =  j  f  d*«p;i(«)e‘^^'“''^.  (517) 

Therefore,  the  spatial  correlation  of  p«i(.R)>  O’**  >* 

C,{RuR7)  =  (p.i(^i)p:i(^j)) 

=  //  ^«// (5-18) 

where 

{pJi(k)P«i*(^'))  =  /(^)/*(<0(f(^*  ~  {K’x~  ktOti^y~~  kyo)),{^<’\^) 

and  /{k)  is  defined: 

/(«)  = 

p2Ku  +  PiK2, 

#2  l2 

l>*o'C3»(l  —  Vq)  +  (1  +  Vo){—kxo>^9  ~  kyoKy  +  — ~  ~  )  •  (5.20) 

Pi  Pi 

Since  f  is  the  Fourier  transform  of  f,  therefore. 


{?(^*  “■  kxOt  ^  ~  ^Vo)f  (^*  ~  kxQ,K^  —  ^yo)) 

”  (2^  /  /  / 1 

g— »'((»*— fc»p)*—(*i—**o)**+(*i~*»p)|f—(*i”'*»o)lfn 


2 

SS  .  *  ■S«(/Cf  ky0fKy  —  kyQ)6^K-^K)f 


(5.21) 
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where  N,(D  is  the  autocorrelation  function  of  and  a]  is  the  standard  deviation 
of  f,  and  S,(k)  is  the  spectnun  of  Nt{f).  In  the  above  derivation,  variable  r  b 
replaced  by  f  =  f  Finally,  substituting  the  above  results  into  Eq.(  5.18),  we 
have 

C.(RuS,)  =  f  j d?ii\fm'S.{K.  (5.22) 

where  .S  =  —  rj  b  the  horbontal  dbtance  between  the  two  receivers;  and  Z  = 

Zx  —  Z]  b  the  vertical  dbtance. 

Assume  Nt(J)  and  its  associated  spectrum  have  one  of  the  following  forms: 


S,(i5) 

S.(«) 


2irll 

-.iUJ-iJil 

g  l«0  ^  J 


(1  —  D  exponential) 


(2  —  U  exponential) 


4lzolyO 

(1  +  (#C*/*o)’)  (1  +  (Wv^vo)*)  ’ 


N,(i)  =  (Gousston) 

5,(k)  = 

we  can  evaluate  £q.(  5.22)by  numerical  integration. 


(5.23) 


(5.24) 


(5.25) 


When  the  two  receivers  in  Eq.(  5.22)  or  in  £q.(  4.9)  of  Chapter  4  are  merged, 
jfi  =  R2,  the  result  b  the  interface  or  volume  scattering  intensity,  respectively. 
Using  these  two  equations,  we  can  examine  the  relative  scattering  intensities  of 
interface  scattering  and  volume  scattering,  given  the  same  average  background 
soxmd  speed  profiles.  The  average  background  profile  for  the  following  examples 
are  two  homogeneous  half-spaces  with  sound  speeds  and  densities  given  in  Fig.5.3. 
Based  on  the  background  profile,  we  can  perturb  the  interface  to  make  it  rough, 
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Figure  5.3:  Background  acoustic  parameters 

and  calculate  the  plane  wave  scattering  intensity  from  roughness  by  £q.(  5.22); 
also,  we  can  keep  the  interface  flat,  but  perturb  the  bottom  sound  speed  in  a  flnite 
layer  immediately  beneath  the  flat  interface,  to  calculate  the  plane  wave  scattering 
intensity  from  volume  inhomogeneities  by  Eq.(  4.9).  By  comparing  the  two  re- 
stilts,  we  will  be  able  to  get  some  quantitative  understanding  of  scattering  from  the 
rough  interface  versxis  scattering  due  to  volume  inhomogeneities.  Fig.5.4  shows  the 
comparison  of  scattered  pressure  magnitudes  vs.  frequency;  two  of  the  four  ciirves 
correspond  to  interface  scattering,  the  other  two  curves  correspond  to  volume  scat¬ 
tering.  In  these  examples,  the  rms  roughness  heights  are  0.1  and  0.2  meters,  and 
the  rms  variations  of  sound  speed  are  5  and  7  per  cent.  In  these  cases,  volume 
scattering  dominates  for  low  frequencies,  and  interface  scattering  dominates  for 
high  frequencies.  It  is  clear  from  these  examples  that  the  roughness  is  a  strong 
scatterer;  however  volume  scattering  cannot  be  ignored  when  the  interface  is  rel¬ 
atively  smooth,  especially  for  low  frequenci^.  We  note  that  the  crossing  points 
between  the  curves  will  move  to  higher  or  lower  frequencies  when  the  roughness 
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height  or  the  volume  randomness  size  changes;  however  the  shapes  of  the  curves 
are  determined  by  the  correlation  functions  of  the  scatterers.  This  suggests  that 
as  long  as  the  roughness  and  the  volume  inhomogeneities  have  different  correlation 
functions,  the  frequency  dependence  of  tne  scattered  field  will  be  different.  There¬ 
fore,  there  is  a  possibility  to  distinguish  whether  volume  scattering  exists  when  we 
know  only  the  roughness  spectra  by  comparing  the  predicted  scattered  field  due  to 
roughness  alone  with  the  total  measured  field. 
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Figure  5.4:  Comparison  of  scattered  wave  magnitudes 


Next,  the  spatial  correlations  of  scattered  fields  will  be  compared.  The  back¬ 
ground  profile  is  the  same  as  in  Fig.5.3.  Eqs.(  5.22)  and  (  4.13)  are  used  to  cal¬ 
culate  the  spatial  correlations  for  rough  interface  and  for  volume  inhomogeneity 
scattering,  respectively.  So  far  the  spatial  correlations  of  scattered  field  are  not 
normalized.  To  examine  the  combined  scattering  effect  of  roughness  and  volume 
inhomogeneity,  we  define  the  normalized  spatial  correlation: 

c(^i,.R;)  =  C{^^,R^)/C{RuR^).  (5.26) 

So  the  correlation  is  normalized  to  unity  at  zero  separation  distance. 

Figs.5.5  and  5.6  show  the  case  when  interface  scattering  dominates.  The  spa¬ 
tial  correlations  in  Fig.5.5  is  not  normalized  in  order  to  compare  the  scattering 
strength;  one  of  the  two  curves  in  Fig.5.5  is  the  spatial  correlation  by  rough  inter¬ 
face;  the  other  curve  is  the  spatial  correlation  by  volume  inhomogeneities.  Notice 
that  the  spatial  correlation  length  for  the  interface  curve  is  about  5  meters,  which 
equals  the  correlation  length  of  the  rough  interface  random  variable  the  spatial 
correlation  length  for  the  volume  curve  is  about  15  meters,  which  equals  the  cor¬ 
relation  length  of  the  volume  random  variable  e.  Fig.5.6  is  the  normalized  spatial 
correlations  of  scattered  field  caused  by  the  roughness,  by  volume  inhomogeneities, 
and  by  the  combined  effect.  Since  the  interface  scattering  is  much  stronger,  the 
curve  for  the  combined  correlation  in  Fig.5.6  is  similar  to  that  of  the  roughness 
correlation.  Therefore,  what  we  would  invert  in  this  case  is  the  roughness  correla¬ 
tion  length. 

Figs. 5.7  and  5.8  show  the  same  type  of  curves  as  in  Figs.5.5  and  5.6  except  that 
neither  scattering  mechanism  dominates.  The  curves  in  Fig.5.8  are  normalized  spa¬ 
tial  correlations  of  the  scattered  field  due  to  interface,  volume  and  combined  effect. 
We  find  from  the  curvature  change  of  the  combined  correlation  that  both  the  5 
meter  correlation  length  from  the  interface  and  the  15  meter  correlation  length 


96 


from  the  volume  can  b .  Identified.  Figs.5.9  and  5.10  are  the  saune  as  Figs.5.7  and 
5.8,  except  that  the  correlation  function  of  the  random  variable  (  of  the  interface 
is  not  Gaussian,  but  exponential.  They  show  similar  resxilts. 

This  suggests  that  by  using  a  series  of  receivers  to  form  spatial  correlations  of  the 
total  scattered  field,  it  may  be  possible  to  establish  whether  the  scattering  process 
is  dominated  by  one  of  the  scattering  mechanisms  or  by  both.  When  one  of  the  two 
scattering  mechanisms  dominates,  we  will  be  able  to  invert  its  correlation  length, 
and  the  contribution  from  the  other  mechzmism  will  be  treated  as  noise;  when 
neither  scattering  mechanism  dominates,  the  correlation  curve  of  the  scattered  will 
have  a  turn  on  its  slope. 
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!  5.5:  Comparison  of  scattered  field  correlations;  interface  roughness  domi- 
without  normalization. 
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Figure  5.6;  Comparison  of  scattered  field  correlations,  interface  roughness  domi¬ 
nates;  normalized. 
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Figure  5.7:  Comparison  of  scattered  field  correlations;  neither  interface  roughness 
nor  volume  dominates;  without  normalization. 
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Figure  5.8:  CompariBon  of  scattered  field  correlations,  neither  interface  roughness 
nor  volume  dominates;  normalised. 
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dominates:  without  normalization. 


Chapter  6 


SPECTRAL 

PARAMETERIZATION  OF 
SCATTERING  FROM  A 
RANDOM  BOTTOM  AND 
POINT  SOURCE 
CONFIGURATION 


The  conventional  parameter  which  quantifies  scattering  from  ocean  bottoms  is 
the  bottom  scattering  cross  section  or  bottom  scattering  coefficient.  Although  this 
parameter  has  been  widely  used,  we  find  it  does  not  have  a  clear  physical  meaning. 
In  this  chapter,  we  examine  the  scattering  concept  and  define  a  new  set  of  param¬ 
eters  to  quantify  bottom  scattering.  In  the  previous  chapters,  we  examined  the 
plane  wave  scattering  problems;  here  we  will  formulate  the  point  source  scattering 
problem  in  terms  of  the  new  parameters. 
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The  physical  concept  of  scattering  is  developed  from  the  analysis  of  wave  prop¬ 
agation  in  a  medium  that  is  homogeneotis  except  for  a  small  region  (the  scattering 
region)  where  the  physical  properties,  e.g.  sotmd  speed  or  density  in  acoustics, 
refractive  index  in  electromagnetic  (EM)  waves,  are  different  from  the  rest  of  the 
medium.  However,  ocean  bottom  scattering  is  different  from  the  classical  scatter¬ 
ing  problem  because: 

(a)  the  averaged  acoustical  parameters  of  the  bottom  are  not  homogeneotis 

and  are  often  layered; 

(b)  the  scattering  region  is  not  confined  to  a  finite  region. 

In  order  to  find  a  set  of  parameters  to  characterize  (parameterize)  this  type  of 
scattering  medium,  it  is  necessary  to  examine  the  classical  definition  of  scattering 
cross  section  which  is  the  parameter  characterizing  the  scattering  region. 

According  to  Morse  and  Feshbach[63]  (see  Fig.6.1),  the  total  wave  function  0  is 
decomposed  into  a  plane  incident  wave  V'«  &  scattered  wave  rj),'. 

V'  =  V'i  +  0..  (6.1) 

The  incident  wave  is  propagating  in  the  positive  z  direction: 

rifi  =  e’*'.  (6.2) 

The  scattered  wave  is  evaluated  at  inffnity,  which  has  the  asymptotic  form: 

/(^,0)(e’*7»');  r-»oo,  (6.3) 

where  r,  and  ^  form  the  spherical  coordinates  with  an  origin  located  in  the 
scattering  region.  The  scattering  cross  section  a  is  defined  as: 
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Figure  6.1:  Geometry  of  a  classical  scattering  problem 

(6.4) 

So,  the  scattering  cross  section  is  a  quantity  .;hat  measures,  at  in&nity,  the  power 
scattered  into  the  direction  per  unit  solid  angle,  per  unit  incident  intensity. 

It  should  be  noticed  that  the  definition  of  or  is  independent  of  source  position 
(plane  wave  incident),  which  is  necessary  for  the  parameterization;  and  a  is  evalu¬ 
ated  at  infinity,  which  means  that  in  the  region  where  a  is  measured,  the  scattered 
wave  hai  no  interaction  'vith  the  scattering  region. 

In  underwater  acousticii,  the  parameters  used  to  characterize  volxime  scattering 
and  interface  scattering  are  volume  and  surface  scattering  coefficients.  They  are 
defined,  respectively,  as 

or.  *  I,[B,  tf>) .  rV(/o(«o.  <h) '  dvy,  (6.5) 

and 

O,  -  /,(«,  •  '■Vl/o(«o,  *)  •  <i5|.  (6.6) 
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Where,  (^o,  ^o)  is  the  incident  plane  wave  direction,  {$,  4)  is  the  scattering  direction, 
/q  b  the  incident  intensity,  /« is  the  scattered  intensity  per  unit  solid  angle  measured 
at  a  distance  r  from  the  scatterer,  dV  and  dS  are  the  insonified  volume  and  area, 
respectively. 

It  is  clear  that  the  definitions  of  and  a,  are  transplanted  directly  from  the 
classical  definition  of  scattering  cross  section  with  the  following  assumptions: 

(a)  the  total  scattering  volume  or  surface  is  divided  into  many  differential  ele¬ 
ments  (small  volumes  or  sxirfaces],  and  each  of  these  differential  elements  is  assumed 
to  be  an  independent  scatterer;  there  are  no  interactions  among  the  scattering  el¬ 
ements.  The  total  scattering  is  the  summation  of  the  contribution  of  each  of  the 
scattering  elements. 

(b)  the  incident  wavefront  should  be  plane  (or  locally  plane  at  the  scattering 
element). 

(c)  compensation  by  multiplying  by  the  factor  r*  to  the  scattering  strength 
means  that  the  scattered  waves  are  assumed  spherical  waves. 

Although  the  scattering  coefficients  have  been  used  extensively,  there  are  con¬ 
ceptual  and  practical  difficulties  when  applying  them  to  ocean  bottom  scattering 
problems. 

First,  the  division  of  the  scattering  volume  or  surface  into  differential  scattering 
elements  is  not  justified  since  interaction  among  the  elements  is  possible.  One  could 
argue  that  for  very  high  acoustic  frequencies,  and  when  the  correlation  length  of 
the  bottom  random  process  is  much  smaller  than  the  linear  dimension  of  each 
scattering  element,  the  interaction  among  the  elements  can  be  ignored.  However, 
a  proper  parameterization  should  not  allow  such  an  assumption. 

Secondly,  the  bottom  medium,  on  average,  is  not  a  homogeneous  medium.  So 
the  incident  wavefront  at  the  scattering  element  may  not  be  planar.  Similarly,  the 
scattered  wave  cannot  be  assumed  spherical,  since  it  can  take  a  curved  path  to  the 
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Figure  6.2:  Multi-path  in  a  layered  structure 

observation  point  via  refraction,  or  divide  into  more  thzui  one  wave  form  via  multi- 
path.  Therefore  the  compensation  for  spreading  loss  by  the  factor  generally  is 
not  justified. 

As  an  example,  Fig.6.2  shows  a  case  where  the  scatterer  is  "sandwiched”  in 
a  sub-bottom  layer.  The  figure  shows  three  possible  incident  paths,  i.e  refracted 
by  the  first  interface,  refracted  first  by  the  first  interface  and  then  reflected  by 
the  second  interface,  and  a  lateral  wave  path.  Also  shown  axe  similar  scattered 
paths.  In  order  to  find  a  natural  set  of  parameters  to  characterize  the  bottom 
scattering  process,  it  is  helpful  to  review  the  parameterization  of  reflection  from 
layered  medium.  Broadly  speaking,  a  layered  medium  itself  is  a  scatterer,  because 
waves  will  be  diffracted  away  from  straight  line  propagation.  Plane  wave  reflection 
coefficients[57]  are  the  parameters  which  fully  describe  layered  media  in  the  sense 
that,  by  knowing  these  parameters,  i.e.  all  the  reflection  coefficients  of  incident 
angles  ranging  from  0  to  90  degrees  and  to  infinite  imaginary  angles  (inhomogeneous 
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Figiire  6.3:  Reflection  from  stratifled  media.  Reproduced  from  Brekhovskih,  Waves 
in  layered  media. 

plane  wave)  [57],  the  influence  of  the  layered  structure  on  wave  behavior  above  the 
layered  structure  is  completely  known.  For  a  particular  sound  sotirce,  the  response 
of  the  layered  medium  is  the  linear  superposition  of  the  reflection  coefficient  at  all 
angles.  For  example,  the  point  source  response  is  given  by  Brekhovskikh  [57](see 
Fig.6.3): 


/(r/a)-»oo 

=  jQ{krsin9)exp\ik{z  +  ZQ)eos$]V  {B)sin9dB .  (6.7) 

It  is  clear  that  in  order  to  obtain  a  set  of  parameters  that  totally  character¬ 
ize  the  random  ocean  bottom,  one  needs  to  consider  the  inhomogeneities  within 
the  bottom,  as  well  as  the  layered  structure  of  the  averaged  background  bottom 
medium  as  one  large  scatterer,  and  to  have  a  plane  wave  impinging  on  to  the  large 
scatterer.  Since  the  scattering  medium  is  random,  the  second  moments  of  the 
randomly  scattered  field  (scattering  strength  and  spatial  correlations  of  scattered 
field)  are  concerned;  hence  the  set  of  parameters  will  be  the  cross  correlation  of  the 
scattering  field  at  the  water/bottom  interface,  induced  by  two  plane  waves  of  unit 
amplitudes  impinging  firom  arbitrary  directions,  including  inhomogeneous  plane 


waves.  This  set  of  parameters  will  be  called  Scattering  Correlation  Coefficients. 
They  form  an  infinite  matrix.  The  diagonal  terms  of  the  matrix  are  correlations 
of  scattering  fields  induced  by  two  plane  waves  impinging  from  the  same  direction; 
the  off  diagonal  terms  are  correlations  induced  by  two  plane  waves  impinging  from 
different  directions.  For  an  arbitrary  sound  source,  the  response  of  the  random 
ocean  bottom  will  be  a  linear  superposition  of  all  the  correlation  coefficients. 

In  order  to  introduce  the  Scattering  Correlation  Coefficients,  the  bottom  is 
modeled  as  shown  in  Fig.6.4,  and  the  single  scattering  approximation  (Extended 
Born  approximation)  is  used.  Although  the  bottom  is  specific,  the  Scattering  Cor¬ 
relation  Coefficients  derived  therefrom  are  general. 

Asstune  that  a  plane  wave  with  horizontal  wave  vector  /Ti  and  unit  amplitude 
is  impinging  on  to  the  bottom: 

=  exp{itrx  •  r  -  inx,z) ,  (6.8) 

where  pxn  is  the  vertical  component  of  the  the  wave  number: 

Hx*  =  yl^x-A-  (6-9) 

The  zeroth  order  solution,  i.e.  the  solution  in  the  absence  of  the  random  quantity 
fc,  is  the  solution  of  waves  in  a  layered  medium,  which  is  available  either  analyt¬ 
ically  or  numerically[4l].  Let  the  zeroth  order  solution  in  the  random  layer  be  a 
complex  quantity: 

r, z)  =  r(Mi, z)e’'^‘  ^  (6.10) 

where  e*'**'^  results  from  Snell’s  law.  This  quantity  is  the  induced  source  strength 
under  the  Bom  approximation. 

In  the  random  layer,  the  wave  equation  governing  the  scattering  field  in  the 
Extended  Born  approximation  is: 
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I  homogeneous 

Figure  6.4:  Schematic  diagram  of  volume  scattering 
(V»  +  kl]p\  S  2A:*£(A)r(/ri,  a)***^* (6.11) 

where  €  =  dc/cj  is  the  random  quantity  due  to  the  random  sound  speed,  is 
the  averaged  wave  number  in  the  layer. 

Since  the  Green’s  fimction  for  a  layered  medium  can  be  written  as[4l]: 

G(.R,  =  ^  /  ^  /  (6.12) 

where  (x,  y,  z)  is  the  receiver  position,  and  (x',  y',  a')  b  the  secondary  source  posi¬ 
tion,  and  G  obeys 

[V’  +  k*]G  =  6{R  -  R%  (6.13) 

the  scattered  field  can  be  expressed  as: 

Pi{R,l!i)  =  (2k\)  j  j^f  €(«)T(iiu^)t>f'^G(R,B;)iv' 
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Changing  the  integration  order,  the  above  expression  becomes: 


(6.15) 


The  scattered  field  in  the  water  column  can  be  written  as  summation  of  up-going 
plane  waves: 


=  £“/ (i’/c-imW.M-Oz"*  V'v'^.  (6.16) 

At  the  water/bottom  interface,  2  =  0,  the  scattered  field  is  continuous;  =  p,. 
Comparing  Eqs.(  6.15)  and  (  6.16),  we  have: 

"(-fi./ri)  =  (7)  //,  /  (6.17) 

Clearly,  m(/cl,/Zl)  measures  the  amount  of  scattering  (including  amplitude  and 
phase)  going  to  /cl  direction  when  a  plane  wave  b  impinging  upon  the  random 
medium  from  /iTi  direction.  Thb  quantity  b  independent  of  source  or  receiver 
position,  and  it  contains  all  the  information  about  the  random  bottom  necessary 
for  determining  the  statbtics  of  the  scattered  field  in  1^  direction. 

To  measure  the  second  moments  of  the  randomly  scattered  field  p\  (the  first 
moment  of  p\  b  zero  because  c  b  assiimed  to  have  zero  average  value),  it  b  enough 
to  know  the  cross  correlations  between  the  m(^,/Il)’s.  The  Scattering  Correlation 
Coefficients  are  defined  as: 


J!>(«i,«2,Mi,;ra)  =  (m(/ri,|ri)  •m*(»Q,/r2)).  (6.18) 

Changing  variables  as  before: 

/  =  -  P*',  /  =  yj[x'  -  x")»  +  (y'  -  iT)’,  (6.19) 

and  considering  the  random  quantity  c  to  have  the  same  statbtics  as  in  the  previous 
chapters: 
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(t(#)<(^”))  =  <r» 


(6.20) 


5(S)  =  J  jN{T^e-**W,  (6.21) 

Kty  fTii  fi2)  can  be  evaluated: 

i?(«l,/ra,/ri,/ra)  =  /<i‘*'ff(«i.*')ff*(«a,a')r(/ri,z')r(Ara,ar') 

=  ZQ{2akiyS{iri  -  Ki)5[(Ari  -  ih)-  (/Ti  -  iTa)] 

J  dz‘g{Kiy3^)g*{K2yz')T{fruz')T*{/r2yz').  (6.22) 

Among  the  four  vectors  in  i?(iq,/Q,  At/Za)  only  three  of  them  are  independent; 
the  fourth  one  b  determined  by  the  Dirac  delta  in  the  above  expression.  This  is  a 
consequence  of  assuming  the  scatterers  are  spatially  homogeneous. 


Knowing  the  Scattering  Correlation  Coefficients,  we  can  construct  the  corre¬ 
lation  of  the  scattered  fields,  induced  by  any  two  incident  plane  waves,  by  using 
Eq.(  6.16): 

{p{{Rufii)p*i{R2,lj^2))  -  J  f  /  /  <^*«2-D(«i,'^,Mi,7ra) 

When  a  specific  acoustic  source  is  used,  the  spatial  correlation  of  the  scattered 
field  can  be  found  by  adding  up  contributions  from  all  incident  plane  waves  with 
appropriate  weighting  functions.  As  an  example,  a  point  source  can  be  expressed 
as  superposition  of  plane  waves[57]: 


(6.24) 


where 
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(6.25) 


VZTrVfcJ-Aij; 

is  the  weighting  function  of  each  plane  incident  wave  for  point  source.  Therefore, 
the  spatial  correlation  of  the  scattered  field  from  a  point  source  is: 

=  //<?«// 

\w(i,u  «)1(p;(A.  wWC*.  /r,)X6.26) 


Combining  Eqs.(  6.23)  and  (  6.26),  we  obtain  the  expression  for  spatial  correlation 
of  the  scattered  field  due  to  a  point  source  as  a  multi-dimensional  integral: 


(p'i(^i,^)pr(^2.<^)> 


j  j  dVi  f  j  <i*/i*2lW^(Mi.^)^*(A*2,2o)] 
j  j  j  j  d^KiDiiCi^K^.iTuirt) 


(6.27) 


When  Ri  ss  Ri,  Eq.(  6.26)  gives  the  scattering  strength;  further,  if  JRi  =  J^, 
it  gives  the  backscattering  strength.  It  is  clear  from  the  general  expression  of 
Eq.(  6.27)  that  all  the  second  moments  of  the  scattered  field  ue  linear  superposi¬ 
tions  of  the  Scattering  Correlation  Coefficients  I?(^i,/^,;ri,/2^). 

There  are  three  points  worth  mentioning.  First,  the  quantity  z')r*(Mj,  s') 
in  Eq.(  6.22)  is  the  incident  strength  at  the  scatterer,  which  generally  is  not  that 
of  a  plane  wave;  it  contains  multipath  and  lateral  wave  contributions. 

Second,  Eq.(  6.26)  shows  that  the  received  scattering  strength  from  a  point 
source  is  the  combined  result  of  scattering  in  all  directions.  It  is  valuable  to  solve 
Eq.(  6.26),  and  to  compare  the  result  with  that  of  the  conventional  definition  of 
scattering  coefficient  in  Eq.(  6.5). 


Realizations  of  the  scattered  field  is  possible  if  we  generate  a  particular  random 
bottom.  The  statistics  of  the  scatterers  can  be  controled  by  specifying  the  spectrum 
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of  the  scatterers.  By  generating  realizations  over  a  nunaber  of  frequencies,  we  will 
be  able  to  synthesize  scattered  time  series  (reverberation  signals) . 
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Chapter  7 


CONCLUSIONS  AND 
RECOMMENDATIONS  FOR 
FUTURE  RESEARCH 


7.1  Conclusions 

Ocean  bottom  scattering  is  a  very  complicated  problem.  This  thesis  is  a  the¬ 
oretical  analysis  of  some  aspects  of  the  scattering  problem.  We  emphasize  on  the 
scattering  from  the  inhomogeneous  volume,  but  the  interface  scattering  is  also  ex¬ 
tensively  covered.  The  four  major  accomplishments  are: 

(1)  Detailed  study  of  spatial  correlations  of  scattered  field.  Spatial  correlations 
of  the  scattered  field  carries  information  which  cannot  be  obtained  from  the  scatter¬ 
ing  intensity.  Conventionally,  bottom  scattering  research  has  been  concentrated  on 
modeling  and  measuring  the  scattered  intensities.  However,  the  scattering  intensity 
only  describes  the  amount  of  energy  being  scattered;  it  cannot  provide  information 
on  how  the  scattered  field  is  spatially  related.  We  found  that  the  spatial  correlation 
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length  of  the  scattered  field  is  related  to  the  statistical  properties  of  the  scatterers. 
By  using  a  wideband  acoustic  source^  and  sweeping  acoustic  frequency,  we  can  in¬ 
vert  the  correlation  length  of  the  scatterers  from  the  measurement  of  the  spatial 
correlation  of  the  scattered  field.  In  addition,  the  infiuence  of  bottom  anisotropy 
is  also  investigated. 

(2)  Introduction  of  the  Scattering  Correlation  Coefficient.  Ocean  bottom  is 
a  complicated  medium,  however  the  scattering  cross  section,  the  major  parame¬ 
ter  quantifying  bottom  scattering,  is  transplanted  from  the  much  simpler  classical 
scattering  problem.  In  many  cases,  especially  in  low  frequency  scattering  problems, 
the  scattering  cross  section  is  not  a  proper  parameter.  By  introducing  the  Scat¬ 
tering  Correlation  Coefficient,  we  have  a  better  understanding  of  scattering  from 
a  random  bottom.  More  importantly,  these  set  of  parameters  has  the  potential 
to  be  used  as  a  standard  quantity  in  comparing  experimental  data  obtained  from 
difierent  sites. 

(3)  Comparison  of  two  diflFerent  scattering  mechanisms.  When  the  water/bottom 
interface  is  not  very  rough,  it  has  been  speculated  that  the  volume  scattering  mech¬ 
anism  might  be  important.  We  combined  the  roughness  and  the  inhomogeneity 
into  one  model  under  the  assumption  that  ihe  roughness  and  the  inhomogeneity 
are  imcorrelated,  and  showed  quantitatively  in  examples  that  the  volume  scattering 
is  important,  especially  for  low  frequencies. 

We  also  showed  that  when  the  scattering  is  caused  by  both  interface  roughness 
and  volume  inhomogeneities,  it  is  possible  to  distinguish  the  contributions  to  the 
scattered  field  from  these  two  scattering  mechanisms  by  examining  the  spatial  cor¬ 
relation  of  the  scattered  field.  To  do  this,  we  have  to  asstime  that  the  roughness 
and  the  volume  inhomogeneities  have  different  statistics;  in  other  words,  they  have 
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to  have  different  spectra. 


(4)  Application  of  the  integro-differential  equation  method.  We  extended  the 
integro-differential  equation  method  to  the  problem  of  weak  scattering  from  a  inho¬ 
mogeneous  layer,  in  which  a  water/bottom  discontinuity  is  included.  By  using  the 
coherent  field  obtained  from  that  integro>differential  equation  as  input  to  calculate 
scattered  field,  we  improved  upon  the  Bom  approximation. 


7.2  Recommendations  for  future  research 

Theoretically,  we  recommend  the  following  research  in  the  future: 

(a)  Develop  a  fast  calculation  technique  to  handle  the  multifold  integration 
problem. 

To  evaluate  the  spatial  correlation  of  the  scattered  field  due  to  a  point  source  is 
of  practical  importance.  Although  we  have  formulated  that  problem  in  this  thesis, 
we  have  not  been  successful  in  solving  it  because  we  do  not  have  an  effective  method 
to  numerically  evaluate  that  multifold  integration  involved  in  that  formulation. 

(b)  Simulate  reverberation  signab. 

Based  on  the  formulation  of  scattering  due  to  a  point  source,  we  can  gener¬ 
ate  Monte  Carlo  type  realizations  of  a  random  bottom  and  calculate  the  scattered 
fields  over  a  band  of  frequencies  for  each  realization.  By  83mthesizing  these  scat¬ 
tered  fields,  we  will  be  able  to  simulate  the  reverberation  signals  in  time  domain 
due  to  a  realistic  finite  duration  pulsing  source.  It  will  help  in  both  designing  ex¬ 
periments  and  processing  real  data. 
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(c)  Application  of  the  Scattering  Correlation  CoeSicient. 

One  goal  is  to  express  the  coefficient  through  measurable  quantities  (ex.  acous¬ 
tic  pressure)  in  order  to  use  it  as  a  standard  to  compare  experimental  data  ob¬ 
tained  from  different  locations.  Another  goal  is  to  use  the  Scattering  Correlation 
Coefficient  to  find  the  applicable  conditions  of  the  classical  scattering  cross  section. 

(d)  Develop  a  model  which  can  predict  scattered  field  when  the  interface  rough¬ 
ness  and  the  volxime  inhomogeneities  are  correlated. 

(e)  Investigate  the  case  when  the  bottom  correlation  function  is  not  a  single 
parameter  form. 

Experimentally,  the  major  issue  Is  what  geological  parameters  are  necessary  to 
be  used  as  inputs  when  we  apply  the  theory  to  experiments.  The  following  are 
some  considerations: 

(a)  The  data  supplied  by  geologists  are  material  properties  and  they  are  not 
directly  suitable  for  acousticians.  Some  of  the  data  are  also  acoustic  properties 
such  as  compressional  velocity  and  attenuation  coefficient;  others  are  indirectly 
related  to  acoustics,  such  as  porosity,  water  content,  grain  weight,  and  grain  size. 
These  materia]  properties  are  especially  important  to  the  acoustic  scattering  pro¬ 
cess.  How  the  variabilities  of  these  parameters  influence  acoustic  scattering  is  not 
clear.  There  has  been  some  research  in  this  area[63],  but  they  are  incomplete.  The 
key  question  is  what  dominates  volume  scattering? 

(b)  In  Chapter  2,  we  modeled  the  scatterers  as  a  variation  of  sound  speed  within 
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the  seabed;  whether  the  variations  of  density  are  also  importzint  as  scatterers  is 
not  clear.  So  we  suggest  that  in  the  geological  surveys  of  an  experimental  site,  the 
correlation  between  sound  speed  and  density  variations  be  conducted.  Ideally,  we 
would  like  to  have  the  joint  distribution  of  the  two  random  variables,  soimd  speed 
and  density.  If  the  contribution  to  scattering  from  density  variations  is  much  less 
than  that  from  sound  speed  variations,  we  can  directly  apply  the  present  theory; 
if  the  density  is  also  proved  to  be  important,  we  will  determine  how  the  density 
variation  is  correlated  with  the  sound  speed  variation.  If  these  two  quantities  are 
perfectly  correlated,  as  Eines[9]  has  assumed,  the  present  theory  can  be  easily  mod¬ 
ified  to  accommodate  the  density  variation  effect;  if  they  are  not  correlated  at  all, 
then  the  scattered  field  would  be  the  superposition  of  the  contributions  from  the 
two  scattering  mechanisms.  In  this  case,  the  theory  in  this  thesis  can  be  improved 
to  account  for  scattering  by  the  two  types  of  independent  scatterers.  However,  be¬ 
cause  there  are  two  random  processes  involved,  pin  pointing  the  contribution  from 
each  mechanism  would  be  difficult.  We  find  from  the  available  data  [53,50]  that 
the  likely  situation  is  that  the  two  random  variables  are  partially  correlated,  but 
there  is  no  quantitative  data  to  infer  the  form  of  the  correlation  function.  Hence 
it  is  of  great  importance  to  have  a  detailed  exploration  of  this  issue. 

(c)  What  is  the  horizontal  variability  of  the  scatterers?  There  exists  a  large 
amount  of  drilling  data  to  infer  the  vertical  variability  of  the  sound  speed  profile 
in  the  bottom,  but  so  far  there  have  been  few  efforts  to  experimentally  determine 
the  horizontal  variability  of  either  soimd  speed  or  density.  Early  models  assume 
that  the  volume  scatterers  are  uncorrelated  |14,46],  but  this  assumption  not  only 
contradicts  intuition,  but  also  disagrees  with  the  drilling  data[50].  The  theory  de¬ 
veloped  in  this  thesis  is  based  on  the  assumption  that  the  sound  speed  variation  is 
well  correlated  horizontally,  and  poorly  correlated  vertically,  because  in  each  sed- 
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imentation  time  period,  similar  material  will  be  deposited  in  one  sediment  layer, 
and  the  sediment  material  will  be  different  for  different  time  periods.  To  get  exper> 
imental  evidence  on  horizontal  variability  of  the  sediment,  we  need  several  groups 
of  three  to  five  core  samples  of  sound  speed  and  density.  Of  primary  importance 
is  the  rough  size  of  the  horizontal  correlation  length.  Since  we  are  interested  in 
low-frequency  scattering  (below  1  kHz),  and  we  can  assume  the  compressional  wave 
speed  is  1500  -  2000  meters/sec.,  the  shortest  acoustic  wavelength  in  the  sediment 
will  be  in  the  range  of  1.5  -  2.0  meters.  If  the  horizontal  correlation  length  is 
shorter  than  one  meter,  we  should  use  the  uncorrelated-scatterer  model;  otherwise, 
the  model  developed  in  the  thesis  should  apply.  One  key  result  in  the  theory  is  that 
the  spatial  correlation  length  of  the  scattered  field  equals  the  horizontal  correlation 
length  of  the  sound  speed  variation;  therefore  knowing  the  true  correlation  length 
of  the  sound  speed  variation  at  the  experiment  site  is  crucial  in  testing  the  theory 
in  an  acoustic  scattering  experiment.  Based  on  the  above  argument,  if  we  choose 
the  highest  frequency  to  be  1  kHz,  and  the  lowest  frequency  to  be  50  Hz,  the  cores 
should  be  separated  from  each  other  by  at  least  2  meters,  and  no  more  than  40 
meters,  and  they  should  not  be  arranged  on  a  straight  line  in  order  to  find  the 
anisotropy. 

(d)  Ideally  we  would  like  to  conduct  the  experiment  in  an  uea  where  the 
seafloor  is  flat.  However,  the  seafloor  will  always  have  some  roughness,  and  the 
roughness  is  an  effective  scatterer  as  shown  in  Chapter  5.  Therefore,  measuring  the 
roughness  spectrum  is  necessary  to  determine  the  contribution  to  scattered  field 
from  roughness.  This  could  be  done  with  joint  registry  of  a  high  frequency  sidescan 
sonar  and  MHz  frequency  narrow  beam  bathymeter. 

A  difficult  and  also  important  item  to  investigate  is  the  geological  link  between 
the  roughness  and  the  variabilities  within  the  volume.  In  Chapter  5,  we  made  a 
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rough  comparison  of  the  scattering  effects  due  to  roughness  and  volume  inhomo¬ 
geneity,  which  is  based  on  the  assumption  that  there  is  no  statistical  correlation 
between  the  two  random  quantities.  If  the  two  quantities  are  related  statistically, 
there  might  be  some  geological  linkage  as  well.  These  kinds  of  information  will 
help  in  improving  the  acoustic  scattering  model.  If  indeed  the  two  quantities  are 
independent,  knowledge  of  the  roughness  spectrum  will  enable  us  to  evaluate  the 
volume  scattering  contribution  by  subtracting  the  roughness  contribution  from  the 
total  scattered  field. 

Another  point  worth  mentioning  concerns  the  dependence  of  the  roughness  spec¬ 
trum  on  the  probing  acoustic  frequency.  The  interface  perceived  by  high-frequency 
acoustic  waves  may  not  necessarily  be  the  same  one  perceived  by  low  frequency 
acotistic  waves.  Therefore,  the  bottom  roughness  spectrum  should  be  measured 
as  a  fimction  of  frequency.  How  to  practically  conduct  measurements  of  this  kind 
remains  to  be  investigated. 

(e)  The  background  acoustic  sound  speed  profile  of  the  seabed  should  be  in¬ 
cluded  in  the  geological  survey,  since  the  amount  of  acoustic  energy  penetrating 
into  the  seabed  heavily  depends  on  the  backgroimd  acoustic  profile,  and  the  amount 
of  acoustic  energy  penetrating  into  the  seabed  is  proportional  to  the  energy  being 
scattered  by  the  volume  inhomogeneities.  The  two  most  important  parameters  are 
the  average  compressional  wave  speed  profile  and  the  average  density  profile.  In 
this  thesis,  we  have  assumed  that  the  average  profiles  of  sound  speed  and  density 
are  constants.  If  they  are  not  constants,  we  should  modify  the  theory  to  accom¬ 
modate  this  effect.  Ideally,  a  bi-static  acoustic  propagation  experiment  should  be 
part  of  a  scattering  experiment  in  order  to  collect  information  on  the  coherent 
component  of  the  field. 
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